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SPATIAL  TRIANGULATION  BI  MEANS  OF  FHOTOGRAMMETHT 


ABSTRACT 


A  method  of  precision  spatial  triangulation  based  on  the 
principles  of  ground  photogramme try  is  outlined.  The  geometry 
and  the  least  squares  adjustment  for  the  orientation  of  an  indi¬ 
vidual  photogrammetric  camera  are  derived  and  the  mathematical 
analysis  for  the  spatial  triangulation  by  jneans  of  intersection 
photo grammetry  is  given.  Numerical  examples  are  added  for  the  major 
steps  of  the  developed  method. 
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I.  INTRODUCTION 


Geodetic  measuring  methods  are  characterized  by  the  measuring  of 
angles  and  distances.  The  corresponding  problems  in  photogramme try  are 
the  determination  of  angles  and  distances  from  photographs.  The  pro¬ 
blem  may  be  considered  in  two  ways,  which  are  basically  different  in 
their  geometrical  approach.  Either  these  quantities  may  be  determined 
analytically  from  the  coordinates  measured,  or  they  may  be  derived  pro¬ 
tectively  by  reproducing  the  bundles  of  rays^  for  each  exposure  of  the 
system  of  cameras.  Only  in  special  cases  does  there  exist  a  simple  re¬ 
lation  between  distances  measured  on  a  photograph  and  on  the  space  ob¬ 
ject  respectively.  In  general,  however,  photogrammetry  deals  only  with 
the  determination  of  angles  and  this  report  will  deal  with  such  measure¬ 
ments  exclusively.  The  general  problem  of  photogrammetric  measurements 
may  therefore,  be  define*!  as  the  triangluation  pf  spatial  positions  from 
angle  measurements.  As  in  geodesy,  one  uses  ip  photogrammetry  either  the 
intersection  of  resectipn  methods. 

The  purpose  of  this  report  is  to  outline  a  method  of  precision 
spatial  triangulation  for  aerial  targets  recorded  on  photographs,  taken 
at  ground  stations.  Therefore,  we  are  dealing  with  a  problem  of  ground 
photogrammetry.  Basic  conditions  of  conventional  ground  photogrammetry, 
as  distinguished  from  aeri&l  photogrammetry,  ares 

1.  The  camera  stations  are  on  the  earth. 

2,  The  camera  stations  do  not  change  their  position  with 
time,  so  that  it  is  justified  to  separate  the  geodetic 
position  measurements  for  the  camera  stations  from  the 
photogrammetric  measuring  procedure. 

However,  one  of  the  most  important  characteristics  of  conventional 
ground  photogrammetry  should  not  be  applied  in  our  case,  i.e.,  the  as¬ 
sumption  that  the  relation  between  the  space  position  of  the  plate  and 
the  local  plumbline  direction,  as  obtained  by  levels,  can  be  used  as  a 
parameter  in  combining  several  stations  for  the  purpose  of  triangulatidh. 
The  physical  reasons  for  this  limitation  are  explained  in  Chapter  II  of 
this  report.  Consequently,  not  all  the  unknown  elements  of  orientation 
can  be  determined  separately  as  is  done  conventionally.  The  elements 
of  the  exterior  orientation  define^, e.  g.,  by  the  azimuth  and  tilt  angles 
of  the  optical  axis  and  by  the  swing  angle  of  the  plate,  must  be  deter¬ 
mined  3 imult ane  oa-aly ,  The  nature  of  the  rigorous  least  squares  solution 
for  the  plate  orientation,  given  in  Chapter  III/3,  makes  it  necessary  to 
determine  these  three  quantities  in  connection  with  the  computation  of 
the  three  elements  of  the  interior  orientation,  denoted  by  the  principal 
distance  and  the  two  plate  coordinates  of  the  principal  point.  These 
complications  call  for  modifying  the  evaluation  procedure  from  the  con¬ 
ventional  ground  methods  t6  those  of  aerial  photogrammetry.  The  latter 
are  characterized  by  the  fact  that  the  determination  of  all  parameters 
of  orientation  must  be  made  simultaneously  because  of  the  changing  posi¬ 
tion  of  the  camera* 
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The  final  triangulation  of  the  recorded  aerial  targets  is  obtained 
from  the  projective  relations  which  exist  between  pairs  and  triplets  of 
photographs  made  from  different  stations.  Thus  the  spatial  triangulation 
is  carried  out  by  intersections  or  resections.  In  the  outlined  method* 
each  station  is  treated  as  an  independent  unit  except  for  the  final  step 
when  corresponding  rays  are  combined  in  a  spatial  triangulation »  The 
present  report  is*  therefore*  limited  to  a  treatment  of  the  problem  based 
on  the  principles  of  "photogrammetry  by  intersection".  Consequently,  the 
photogrammetric  measuring  procedure  is  applied  at  each  measuring  station 
independently.  The  results  for  each  recorded  target  point  are  expressed 
in  parameters  which  are  related  only  to  the  particular  measuring  station, 
e»  g.jby  two  position  angles. 

This  method  is  not  the  only  photogrammetric  method  for  triangulating 
the  spatial  position  of  a  point.  As  a  matter  of  fact,  a  more  elegant  and 
more  economic  method  may  be  used.  This  is  based  on  the  projective  rela¬ 
tions  between  photographs  taken  from  different  locations,  which  may  not 
only  serve  to  triangulate  recorded  points,  but  present  the  possibility  to 
determine  the  relative  orientation  of  these  photographs  exclusively  from 
their  projective  properties.  In  connection  with  such  a  method,  it  is 
possible  to  measure  two  plates  simultaneously  with  stereoscopic  reading 
devices.  These  methods  will  lead  to  a  mechanical-optical  solution  of  our 
problem  with  the  aid  of  a  high-precision  stereo  plotting  machine.  This 
is  essentially  a  three-dimensional  stereoscopic  comparator  by  which  the 
spaLtial  triangulation  problem  is  simulated. 

Although  this  paper  deals  only  with  the  method  of  photogramraetjy  by 
intersection,  the  mathematical  analysis  of  the  problem  of  plate  orienta¬ 
tion  is  presented  in  a  form  which  may  be  suitable  for  future  stereophoto- 
grammetric  measuring  methods. 


II .  PHOTOGRAMMETRIC  BY  INTERSECTION 

The  spatial  triangulation  of  a  recorded  point  from  photographs  taken 
at  two  ground  stations  may  be  computed  in  ground  photogrammetry  by  the 
intersection  method.  (Fig.  1)  The  horizontal  directions  and  the  eleva¬ 
tion  angles,  necessary  for  the  intersection,  are  obtained  if  th£  princi¬ 
pal  distances  d  are  known,  by  measuring  the  plate  coordinates  (x  and  y) . 

The  directions  of  the  camera  axes  against  the  base  line  or  any  other  refer¬ 
ence  datum  are  obtained  either  directly  by  readings  of  horizontal  and 
vertical  circles,  or  indirectly  with  the  aid  of  additional  reference  points, 
whose  space  coordinates  must  be  known.  The  use  of  dials  so  far  has  not 
proved  ho  be  sufficiently  accurate  to  eliminate  the  need  for  control  points. 
Besides  the  unavoidable  systematic  instrumental  errors,  such  (l  method  is 
affected  by  deflections  of  the  vertical  at  the  stations,  since  the  setting 
of  the  phototheodolit©  (photogrammetric  camera  plus  theodolite)  is  obtained 
with  reference  to  the  vertical  as  indicated  by  levels.  Precision  photo- 
grammetry  should, therefore, use  the  data  of  exterior  orientation  obtained 
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from  dial  readings  with  or  without  correction  for  plumbline  anomalies 
as  approximation  values  only  and  evaluate  the  final  orientation  of  the 
plate  with  the  aid  of  recorded  control  points.  It  is  evident  that  the 
coordinates  of  such  reference  points  must  be  recorded  in  a  way  such 
that  they  can  be  accurately  identified  and  measured.  It  should  be  men¬ 
tioned  here  that  the  use  of  control  points  in  geodetic  practice  is  con¬ 
sidered  essentials  especially  in  order  to  eliminate  the  influence  of 
systematic  instrumental  errors  on  the  final  plate  orientation.  The  fact 
that  it  is  possible  to  compute  from  recorded  control  points  corrections 
to  the  orientation  elements  of  the  plate*  shows  that  the  elements  of 
orientation  may  be  obtained  exclusively  from  measurements  of  such  con¬ 
trol  points.  Hence*  dials  and  levels  are  not  basic  elements  of  photo- 
grammetric  instrumentation  design.  They  may  be  useful  to  increase  the 
computing  economy  but  it  should  be  understood  that  the  parameters  ob¬ 
tained  from  such  auxiliaries  should  be  considered  only  as  approximations. 
Consequently*  a  theory  for  a  photogrammetric  measuring  method  for  pre¬ 
cision  triangulation  must  be  established  for  a  completely  independent 
photogrammetric  camera.  The  so-called  "Ballistic  Camera" 9  which  is  used.) 
e.  g.^to  record  ballistic  data  of  full-scale  missiles  in  flight*  repre¬ 
sents  this  type  of  photogrammetric  camera. 

Summarizing*  the  interior  and  exterior  orientation  of  a  camera  may 
be  computed  in  a  given  coordinate  system  from  recorded  control  points 
if  the  coordinates  of  such  points  are  given  within  the  same  reference 
datum.  Hence*  the  direction  to  an  additional  recorded  target  point  may 
be  obtained  from  the  computed  plate  orientation  and  the  measured  plate 
coordinates  of  the  point.  The  result  may  be  expressed  by  two  position 
angles* e.  g.*by  an  azimuth  and  elevation  angle.  If  such  results  from 
two  or  more  stations  are  combined*  the  space  position  of  the  target  point 
may  be  found  by  intersection. 

The  determination  of  space  positions  of  missiles  from  photographs 
taken  at  several  ground  stations  may  be  considered  as  a  similar  problem. 
However*  this  problem  is  different  in  that  the  camera  axes  are  directed 
into  space  so  that  in  general  no  terrestrial  reference  points  will  be 
within  the  angle  of  view  of  the  camera  lens  or  at  least  they  will  be 
recorded  only  on  the  edge  of  the  plate.  In  such  a  case,  it  is  possible 
to  use  stars  as  control  points.  However*  this  procedure  will  give  rise 
to  some  difficulty  which  will  be  considered  next.  He  have  seen  that  the 
control  points  must  be  known  in  the  same  coordinate  system  in  which  the 
final  orientation  of  the  plate  and  the  final  measuring  data  are  required. 
This  makes  it  necessary  to  establish  the  relation  between  the  right  ascen¬ 
sion  and  declination  coordinates  of  the  stars  and  the  local  earth  fixed 
coordinate  system  of  the  measuring  stations.  It  is  well  known  that  there 
are  formula#  which  allow  us  to  express  the  azimuth  and  elevation  angles 
of  a  star  as  functions  of  right  ascension*  declination*  sidereal  time 
and  geographic  latitude  of  the  station.  The  right  ascension.;  and  declina¬ 
tion  values  as  well  as  the  time  measurements  may  be  obtained. (With  suf¬ 
ficient  accuracy.  Hence*  the  problem  is  to  determine  the  geographic 
coordinates  of  the  stations.  The  determination  of  the  latitude  and  longi¬ 
tude  of  the  stations  may  be  obtained  by  astronomical  observations .  These 
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results,  however,  are  not  useful,  due  to  the  fact  that  the  measurements 
are  affected  by  vertical  deflections,  caused  mainly  by  mass  irregularities 
in  the  crust  of  the  earth.  The  same  problem  exists  in  the  astronomical 
determination  of  the  azimuth  of  a  baseline.  Another  source  for  obtaining 
the  geographic  coordinates  are  the  triangulation  results.  However,  tri¬ 
angulation  systems  must  be  considered  as  local  constructions  distorted 
by  unknown  amounts,  Especially  poor  is  the  | orientation  of  extended  tri¬ 
angulations.  In  addition  the  assumption  is  not  justified  that  the  primary 
orientation  of  the  reference  ellipsoid  is  so  good  that  the  relation  be¬ 
tween  the  celestial  and  the  earth  fixed  system  can  be  obtained  accurately 
enough  to  serve  as  a  basic  parameter  for  a  precision  triangulation  by  two 
or  more  photogrammetric  stations. 

In  the  following,  a  method  is  developed  in  which  this  difficulty 
will  be  eliminated.  It  is  obvious  that  the  orientation  of  the  plate  may 
be  obtained  in  the  right  ascension  -  declination  system  from  star  measure¬ 
ments.  With  the  measured  plate  coordinates  for  any  recorded  aerial  target 
point,  the  position  of  such  a  point  in  the  sky  may  be  expressed  by  right 
ascension  and  declination  angles,  thus  treating  the  recorded  point  as 
an  artificial  star.  Since  the  stars  are  at  essentially  infinite  distance, 
a  given  portion  of  the  sky  will  show  an  identical  pattern  independent  of 
thb  location  of  the  observer.  However,  the  apparent  place  of  the  target 
point  in  the  sky  will  change  with  the  position  of  the  observer.  There¬ 
fore,  the  target  point  will  have  a  different  set  of  right  ascension  and 
declination  values  for  each  measuring  station.  Combining  the  results  of 
two  stations,  it  is  possible  to  determine  the  spatial  angle  between  two 
such  directions  (angle  a  in  Fig.  2).  If  a  ttiird  measuring  station  is 
combined  with  each  of  the  two  other  stations,  three  such  spatial  angles 
will  be  determined.  By  combining  the  angles  with  the  slant  distances 
between  the  measuring  stations,  which  may  be  computed  from  local  geodetic 
data,  a  pyramid  is  formed.  We  now  have  a  typical  problem  of  photogram- 
metry,  namely,  that  of  determining  the  coordinates  of  a  point  in  space 
from  three  pdsition  angles  with  reference  to  three  known  points,  i.e., 
the  well-known  problem  of  resection  in  space.  With  this  step  the  tri¬ 
angulation  problem  is  solved  without  resorting  to  the  critical  values 
of  geographic  coordinates.  The  result  is  obtained  by  making  the  transi¬ 
tion  from  the  celestial  to  the  earth-fixed  system  at  that  stage  in  the 
computing  procedure  when  a  quantity  is  obtained  ..whos  e  magnitude  is  inde¬ 
pendent  of  the  coordinate  system  used  ,  The  spatial  angles  between  the 
determining  rays  of  the  target  point  are  such  independent  parameters. 

It  is  noted,  that  if  the  exposures  are  made  synchronously  at  all  measuring 
stations,  the  time  does  not  enter  into  the  solution,  with  the  exception 
that  an  approximate  sidereal  time  is  needed  to  compute  the.  refraction. 
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Figure  2 


If  the  star  exposures  at  the  different  stations  are  made  at  different 
times,  this  time  difference  only  must  be  known  accurately „  This  measure¬ 
ment  will  enter  into  the  solution  as  a  corresponding  correction  to  the 
right  ascension  values. 

The  conclusion  is  that  for  the  establishment  of  a  precision  photo- 
grammetric  measuring  system,  using  stars  as  control  points,  three  stations 
are  sufficient  and  necessary  in  ortfer  to  obtains  rigorous  solution  for  a 
spatial  triangulation.  A  triangulation  by  intersection  from  any  two  sta¬ 
tions  may  be  considered  as  adequate  only  if  a  systematic  orientation  error 
of  the  base-triangle  has  been  determined  from  the  results  ofa  former 
rigorous  three-station  Solution.  Thus,  a  correction  to  the  azimuth  of 
the  local  triangulation  system  would  be  determined.  This  procedure  will 
be  referred  to  in  this  report  as  .the  ®calibration  of  a  photogrammetric 
base  line®. 


III.  ORIENTATION  OF  A  PHOTOGRAMMETRIC  CAMERA 


1.  The  photogrammetric  camera  and  the  perspective  properties  of 
the  photograph. 

In  a  previous  chapter,  it  has  been  shown  that  a  precision  photo¬ 
grammetric  triangulation  method  should  be  based  on  measurements  obtained 
from  an  independent  photogrammetric  camera,  i.e.,  measurements  not  ' 
dependent  upon  azimuth  and  elevation  dial  readings.  Mathematically, 
we  have  to  consider  only  the  optical  system  as  Hie  center  of  projection 
and  the  photographic  .plate  as  part  of  a  plane  cutting  a  bundle  of  per¬ 
spective  rays.  Each  perspective  may  be  explained  as  a  figure  obtained,  ., 
by  cutting  a  bundle  of  perspective  rays  with  a  plane.  Consequently, 
each  photograph  represents  a  perspective  figure.  We  distinguish  between 
the  dispositive  if  the  plane  is  situated  between  the  object  points  and 
Idle  center  of  projection,  and  the  negative,  if  the  center  of  projection 
is  located  between  the  object  points  and  the  image  plane. 

The  relative  situation  of  the  individual  perspective  rays  of  such 
a  bundle  is  uniquely  determined,  and  it  may  be  congruently  reconstructed 
if  the  distance  of  the  perspective  center  (in  practice,  the  rear  nodal 
point)  from  the  image  plane,  i.  e.,  the  principal  distance  d,  and  the 
location  of  the  principal  point  on  the  image  plane  are  known.  The  H prin¬ 
cipal  point®  is  the  intersection  of  a  line  perpendicular  to  the  plate 
through  the  perspective  center.  This  line  is  called  the  photograph  per¬ 
pendicular.  The  location  of  the  principal  point  on  the  plate  may  be 
given  by  two  rectangular  coordinates  in  any  plate  coordinate  system  repre¬ 
sented  by  certain  fiducial  marks,  and  be  denoted  by  A  x  and  a  y,  (Fig.  3) 
The  principal  distance  and  the/ coordinates  of  the ' principal  point  are  the 
three  elements  of  interior  orientation.  A  photograph  of  which  the  interior 
orientation  is  known  is  called  a  photogram. 
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Figure  3 


Thus,  the  position  of  the  center  of  projection  is  determined  by  three 
perspective  rays.  The  same  three  rays  determine  the  reciprocal  position 
of  the  center  of  projection  with  respect  to  the  corresponding  point-objects. 
The  three  space  coordinates  of  the  center  of  projection,  (X^.*  Yj,  Z^).  the 
spatial  direction  of  the  plate  perpendicular  (optical  axis)  expressed  e.g.^ 
by  two  position  angles  (A  and  y)  and  the  swing  angle  CO  of  any  reference  / 
line  on  the  plate,  represent  the  elements  of  exterior  orientation.  (Fig.  3) 
Hence,  the  exterior  orientation  is  determined  by  six  parameters. 

2.  Mathematical  Analysis  of  the  Orientation  Problem 

In  the  preceding  chapter,  we  have  seen  that  the  interior  orientation 
of  a  camera  calls  for  the  determination  of  three  parameters  and  that  the 
exterior  orientation  requires  six  elements.  In  our  case  the  space  posi¬ 
tion  of  the  center  of  projection  may  be  determined  separately  from  geodetic 
measurements,  thus  reducing  the  number  of  unknown  elements  for  the  ex¬ 
terior  orientation  to  three.  Hence, the  absolute  orientation  of  a  plate 
in  our  case  calls  for  the  determination  of  six  elements,  three  elements 
of  interior  and  another  three  elements  of  exterior  orientation.  These  six 
parameters  of  the  absolute  orientation: \ 

(1)  Principal  distance  -  d 

(2)  x-coordinate  of  the  principal 
point  -  Ax 

(3)  y*-coordinate  of  the  principal 
point  -  Ay  , 

(U)  the  a2imuth  angle  of  the  plate 
perpendicular  -  A 

(5)  the  tilt  angle  of  the  plate 
perpendicular  -  v 

(6)  the  swing  angle  of  the  plate 
coordinate  system  -  K 


The  three  elements 
of  interior 
orientation 


The  three  elements 
of  exterior 
orientation 


We  have  seen  that  the  orientation  of  the  plate  must  be  determined 
from  recorded  stars  as  control  points.  Therefore,  the  elements  of  ex¬ 
terior  orientation  are  determined  with  reference  to  the  right  ascension- 
declination  system. 

The  right  ascension  and  declination  coordinates  are  rectangular 
spherical  coordinates  and  may  be  represented  on' a  unit  sphere.  In  order 
to  relate  these  spherical  reference  coordinates  to  the  plane  plate  coordi¬ 
nates,  a  transformation  of  one  of  the  two  systems  is  necessary.  With 
regard  to  future  theoretical  work,  it  was  decided  to  transform  the  spheri¬ 
cal  star  coordinates  into  plane  coordinates  by  projecting  the  stars  in  a 
plane, tangent  to  the  unit  sphere.  The  plane  coordinates  of  the  stars  are 
obtained  according  to  the  principle  of  central  projection  with  the  point 
of  projection  in  the  center  of  the  sphere  (Fig.  U),  If  the  plane  is  tangent 
at  the  celestial  pole  and  the  coordinate  system  is  oriented  in  such  a  way 
that  the  £ -axis  represents  the  celestial  meridian  through  the  origin  of 
the  right  ascension  measurements  and  the  h  -axis  is  perpendicular  to  the 
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Figure  U 


4;  -axis  3  the  plane  (standard)  coordinates  £  and  ^  for  a  star  may  be  ex¬ 
pressed  as  functions  of  right  ascension  (RA)  and  declination  (6).  The 
formulas  may  be  read  directly  from  Fig.  li.  They  ares 


+  £  -  -cot  6  cos  RA 
r  r 

*h  -  +cot  6  sin  RA 
l  *  r  r 


(1) 


The  subscript  "r"  indicates  that  the  value  has  been  corrected  for  refrac¬ 
tion.  The  astronomical  refraction  for  the  stars  may  be  computed  as  a 
function  of  zenith  distance  with  the  conventional  formulas  used  in 
astronomical  and  geodetic  practice .  The  effect  of  refraction  on  the 
right  ascension  and  declination  of  a  star  may  be  obtained  with  well- 
known  formulas  derived,  e.  g., in  "Textbook  on  Spherical  Astronomy"  by 
W.  M.  Smart,  or  in  "Elements  of  Practical  Astronomy"  by  W.  W.  Campbell. 


Depending  on  the  geographical  position  of  the  camera  station,  it 
may  sometimes  become  necessary  to  choose  a  plane  tangent  to  the  celes¬ 
tial  equator,  preferably  at  the  point  of  origin  of  right  ascension.  In 
addition,  we  will  have  need  for  a  plane  tangent  at  the  zenith  of  the 
camera  station.  In  such  a  case,  the  £  -axis  represents  the  meridian 
through  the  measuring  station.  The  ^  -axis  is  perpendicular  to  the-^  - 
axis.  The  standard  coordinates  are  in  such  a, cases 


♦  £  (north) 


+  ^  (east) 


tan  zr  cos  A 


tan  z  sin  A 
r 


tan  6  -  tan  f£  cos  t 

r  r 

br 


_ r  "  '  r 

cos  t  +  tan  &  taii  5 

T*  *  T 


(2) 


-  sec  i  sin  t 

r _ r _ 

cos  t  +  tan  &  tan  6 
r  r  r 


z  =  Zenith  distance  of  the  star 

.  A  =  Azimuth  of  the  star,  counted 
clockwise  from  the  south. 

$  =  Geographic  latitude  of  the 
camera  station 

t  »  Hour  angle  -  Sidereal  time  of 
exposure  (9)  minus 
Right  ascension 
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We  have  seen  that  a  photograph  may  be  considered  as  an  exact  central 
projection,  (The  necessary  corrections  to  the  measured  plate  coordi¬ 
nates  due  to  distortion  will  be  discussed  later).  Therefore,  a  star 
photograph  represents  a  central  projection  of  a  certain  portion, of  the 
sky..  Similarly,  the  projection  of  the  stars  in  a  plane  tangent  to  a  unit 
sphere  is  an  exact  central  projection.  Both  projections  may  be  assumed 
to  have  the  point  of  projection  in  common.  Hence,  the  plate  images  de¬ 
noted  by  the  plate  coordinates  x  and  y  and  the  corresponding  star  pro¬ 
jections  denoted  by  the  standard  coordinates  ^  and?^  ,  may  be  considered 
as  photographs  taken  from  the  same  point.  We  must  now  determine  the 
geometric  relations  which  exist  between  two  photographs  taken  from  the 
same  point.  In  accordance  with  the  two  different  geometric  approaches 
to  the  problem,  there  are  two  different  mathematical  analyses. 

First  we  will  consider  the  projection  method,  which  is  felt  to  be 
more  intuitive,  because  the  unknowns  of  the  solution  have  precisely 
defined  physical- meanings. 

Solution  a) 


In  Fig.  (£),  L  is  the  center  of  projection.  The  problem  is  to 
orient  the  plane  containing  the  plate  which  is  shown  as  a  diapositive, 
in  such  a  way  that  a  bundle  of  rays  originating,, from  L  and  passing 
through  the  plate  images  (S)  intersects  a  plane  tangent  to  the  unit 
sphere  at  the  point  L*  in  the  corresponding  star  projections l,  (S1), 

We  introduce  the  following  coordinate  systems g 


Standard  coordinates.  Plane  rectangular  coordinates  of 
a  star  in  the  plane  of  p£  oj  ebl&om^  ofljih at*  .the 

point  L‘,  in  which  the  projection  plane  is  tangent  to 
the  unit  sphere.  ^  represents  the  meridian  on  the  unit 
sphere  through  the  point  of  origin  of  the  azimuth  measure¬ 
ments  . 


x,  y 


Primary  plate  coordinates.  Plane,  rectangular  coordinates 
of  an  image  point,  related  to  a  system  whose  origin  (0)  is 
the  intersection  of  the  x  and  y  axes,  established  by  fidu¬ 
cial  marks . 


x,  y 


x',  y' 


Oriented  plate  coordinates.  Plane,  rectangular  coordi¬ 
nates  of  an  image  point.  The  y-axis  is  the  line  of  inter¬ 
section  of  the  image  plane  with  a  plane  perpendicular  to 
the  plane  of  projection  which  contains  the  plate  perpen¬ 
dicular.  The  origin  is  the  isocenter  X.  The  x,  y  system 
is  tilted  against  the  x,  ,y  system  by  the  swing  angle 
( counted  counter-clockwise) 


Plane  rectangular  coordinates  of  a  star  in  the  plane 
of  projection.  The  systejn' corresponds  in  its’  orientation 
to  the  x,  y  system.  The  origin  is  the  corresponding 
isocenter  . 


1 6 


f>/one  '■  off>/~ojec//o*7 


Figure  5 
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4  x  and^  y  denote  the  coordinates  of  the  principal  point 
(P)  in  the  x,  y  system. 


A  »  azimuth  of  the  plate  perpendicular 
v  *  tilt  of  the  plate  perpendicular' 
d  *  principal  distance 

x  =  (x  »^x)  cos (y  ~^y)  sin^ 
y  =  (y  -^y)  cosAf  +  (x  -Ax.)  sin K  +  d  •  tan  ^ 


From  Figures  £  and  6S  we  obtain 


HI  =  HL  - 


sin  v 


*  h 

01  -  OL 


(U) 


sin  v 


A GLS  ****AliS  L.  therefore 
o  o 


0  S. 


o  ® 


QL 


HL 

HS~ 


(5) 


and 


x_ 

x 


LS 


rr 


GS  ' 
o 

■nr 


GL 

H§~ 


(6) 


I  1  ,  ,  7 

Since  HS  “  LH~y  and  GS  ■  GL+y  ,  we  have,  from  formulae  U,  5,  6, 


GLJy  hy  a  .  '  m  QL°x  hx 

ffiT-y  d-y  sin  v  HL-y  d-y  sin  v 


(7) 
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Furthermore % 


^  ®  -  x  sin  A  -  (  y  +  h  tan  ^  )  cos  A 

.  i  ».  i 

V  ■  -(y  +  h  tan  •£•)  sin  A  +  x  cos  A 
For  our  case  h  =  1  and  therefore 


*  d-y  sin  v 


and 


d-y  sin  v 


>  .%  V  \  * 

5  °  -  (y  +  tan  ^  )  cos  A  -  x  sin  A 

^  =  ■»  (y  +  tan  sin  A  +  x  cos  A 

From  (7  )  and  (8  )  we  derive 


4 


+y  cos  A  +  x  sin  A 
-  d  +  y  sin  v 


y  sin  A  -  x  cos  A 
-  d  +  y  sin  v 


-  tan  cos  A 


-  tan  s-  sin  A 


(8) 


(7’) 


(8*) 


(9) 


Substituting  (3)  into  (9) 

‘{("(y-'Ay)  cos/C  +  (x-^x)sin&D  cos  v+d  sin  v}  cos  A+f(x-Ax)cos/f -(y-idy)  gin/3  sin  A 


+  [(y-  A  y )  cos  K  +  (x-4  *)  sinjJfj  sin  v  cos  v 


(I0a) 


^  COB*+  (^)3^<jcos  v+d  sin  v^sin  A-^(x~Ax)cos<-(y-^y)aiiV^cos  A 

+£(y-Ay)cos<+(x-iAx)sin*f]  sin  v-d  cos  v 
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The  formulas  (lOa)  are  in  agreement  gith  the  wellknown  formulas 
Nos*  12a  and  12b  derived  by  0.  v,  Gruber, 

The  formulas  (lOa)  express  the  standard  coordinates  4  and  in 
terms  of  the  measured  plate  coordinates  x  and  y  in  any  rectangular 
plate  coordinate  system  and  in  the  six  elements  of  orientation,  namely, 
the  three  elements  of  interior  orientation  d,dx,  and*4y  and  the  three 
elements  of  exterior  orientation.  A,  v,  and  %  .  These  six  parameters 
are  the  unknowns  of  the  solution.  Because  each  star  gives  rise  to  two 
equations,  one  for  £  and  one  for  7^  ,  three  stars  are  necessary  and  suf¬ 
ficient  to  obtain  a  unique  solution.  This  result  is  in  agreement  with 
the  previously  stated  principle  that  the  center  of  projection  is  fixed 
with  respect  to  the  plate  as  well  as  object-points  by  any  three  perspec 
tive  rays. 

Solution  b) 


In  the  preceding  paragraph,  the  spatial  direction  of  the  optical 
axis  within  the  chosen  reference  system  was  expressed  by  two  position 
angles,  namely  an  azimuth  angle  (A)  and  a  tilt  angle  ( v )  (Fig.  $) . 
Although  these  parameters  are  essential  for  the  final  triangulation 
problem,  such  an  interpretation  of  the  problem  sometimes  causes  diffi¬ 
culties  due  to  the  fact  that  the  concept  of  an  azimuth  angle  does  not 
exist  anymore  if  the  tilt  angle  v  becomes  zero.  In  order  to  have 


Ferienkurs  in  Photogramme trie  by  0.  v  =  Gruber,  Stuttgart-Verlag  von 
Konrad  Wittwer,  1930. 
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a,  ^-tilt,  denotes  the  angle  between  the Plan®  and  a  plane 
parallel  to  the  Zf  -axis  which  contains  the  optical  axis,  to, ^ -tilt, 
denotes  the  angle  formed  by  the  optical  axis  and  its  projection  in  the 
||  -plane.  The  relations  between  the  a,co,  Wangles  and  the  A,v  and 
\K)  angles  respectively,  are; 


sin  co  ■  -  sin  v  sin  A 
tan  a  =  -  tan  v  cos  A 
cot =  cos  v  tan  A 


cos  v  »  cos  a  cos  co 
tan  A  »  tan  co  cosec  a 
K  ) 

where  (tf)  ■=  swing  angle  of 
method  a. 


The  relations  between  the  standard  coordinates  4  and  and  the 
measured  plate  coordinates  x,  y  may  be  obtained  by  applying  the 
formulas  (5)  twice.  The  result,  obtained  by  a  double  projection,  is 
as  follows? 


|“(y-^y)siny»- (x-^x)oDSjjcosa+|dcos  <o-^(y-Ay)cos<-(x-Ax)sin*|sin  cousin 
■|4  cos  co-  (y-Ay)  cos<-(x~4x)  sin&Jsin  cojcos  a- £(x-4x)  cosK+ (y-£y)  sin*] 


sm  a 


d  sin  co* £(y-Ay)  cos/6*(x-Ax),®in  /fj 

COS  <0 

fd  cos  co—  I 

Ry-4y)  cos<-(x-4x)  sinxjsin.coj 

jcos  a-j 

f(x-Ax)  cos4f+  (y-4y)  sinX 

[sin  a 

These  formulas  agree  with  Gruber's  formulas  Nos.  13a  and  13b 
Solution  c) 


*. 


We  will  now  consider  the  mathematical  analysis  of  the  orientation 
problem  based  on  the  principles  of  photogrammetfy  by  analytical  means. 

In  Fig.  (8),  L  is  the  center  of  projection.  A  bundle  of  rays  is  cut  by 
two  planes.  One  plane  contains  the  plate  (negative),  th^  other  plane 
represents  a  plane  tangent  to  the  unit  sphere  at  point  L1  with  the  center 
at  L.  A  star  image  in  the  plane  containing  the  plate  i^  denoted  by  S. 

The  corresponding  image  in  the  plane  of  projection  is  S  .  We  again  intro¬ 
duce  the  following  coordinate  systems? 


4  >  n 


Standard  coordinates.  Plane  rectangular  coordinates  of  a 
star  in  the  plane  projection.  The  origin  is  at  the  point 
L'j  on  which  the  plane  is  tangent  to  the  unit  sphere.  £ 
represents  the  meridian  on  the  unit  spheri  through  the  7 
point  of  origin  of  the  azimuth  measurements. 


*  See  reference  note  on  preceding  page. 
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Figure  8 


x,  y  Primary  plate  coordinates.  Plane  rectangular  coordinates 
of  an  image  point,  related  to  a  system,  whose  origin  (0) 
is  the  intersection  of  the  x  and  y  axes,  established  by 
fiducial  marks. 


x,y  Oriented  plate  coordinates.  Plane, rectangular  coordinates 
of  an  image  point.  The  y4axis  is  the  line  of  intersection 
of  the  image  plane  with  a  plane  perpendicular  to  the  plane 
of  projection  which  contains  the  plate  perpendicular.  The 
originjLs  the  isocenter  I.  The  x,  y-system  is  tilted  against 
the  x,y»system  by  the  swing-angle K (^counted  clockwise). 

x’ ,y’  Plane rectangular, coordinates  of  a  star  in  the  plane  of 
projection.  The  system  corresponds  in  its  orientation 
to  the  X,y  system.  The  origin  is  the  corresponding  isocenter 
I* .  x  and  y  denote  the  coordinates  of  the  origin  of  the 
x,y  system  (0)  in  the  x,y  system. 


A  =  azimuth  of  the  plate  perpendicular 
v  ■  tilt  of  the  plate  perpendicular 
d  =  principal  distance. 


From  Figures  8  and  9,  we  obtain  again s 


HI  ■  HL  ■ 


sin  v 


d» 


01 


OL 


h _ 

sin  v 


AQLS  AHS  L,  therefore 


“o' 

“3E“ 


HL 

H5 

o 


and 


i 


“o 

a  mr~ m 


QL 

^o 


Since  HS  -  LH-y  and  GS1 
o  o 


GL+y  we  obtain  again,  from  formulas  It,  £  and  6, 


7 


hy 

d-y  sin  v 


and  x 


hx 

d-y  sin  v 


2U 


T6e 
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Furthermore t 

.1  i  . 

|  -  x  sin  A  -  y  cos  A 

it 

7^  »  -  y  sin  A  +  x  cos  A  + 

where  are  the  standard  coordinates  of  I*  in  the ^ system,  and 

x  =  X  cos  )C  +  y  siny^*  x 


(12) 


y  =  y  cos  /C  -  x  sin/{  +  yQ 

where  x  ,  y  are  the  coordinates  of  0  in  the  x,y  system* 

For  h  =°1  a8d  by  means  of  formula  Tq  we  obtain  from  formula  11s 


i  .  ».i.a  *  :  y  g»  *  *  yj. 

»  d  -  y  sin  v  '■» 


(13) 


7\  = 

-  y  sin  A  +  x  cos  A 

d  -  y  sin  v 

and  from 

Formula  (l2)s 

a^x  +  b^y  +  c^ 

•a  " 

a  x  +  b  y  +  1 
o  o^ 

- 

agx  +  bgy  +  cg 

i_  i  _ 

a  *  +  b  y  +  1 
o  o  * 

where 

t 

-h 

sin  (A-4J)  +4  3in  & 

®1  “ 

t 

d  -  y_ 

(1U) 


-h  cos  (A-4$  -4  cos  A* 

t 

d  -  y 

*  A 


+h  cos  (A-Ad  +)^  sin.jt 


-  h  sin  (A-K)  cos  ifc 
d'  - 
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-h'(x  sin  A+y  cos  A) 

o  •'o  * 

°1 - dJ-TT- - 


sin  K 

v-d<- 

b  .  .  2°l£_ 

o  d'  -y 


-h* (y  sin  A»x  cos  A) 
wo _ o  _ 

df  -  y 

*  n 


a  ..  . 

o  d*-y„ 


b  .=  .  glJL 

o  d1  -y„ 


0  -  -  ■*  o 

and  the  elements  of  orientation  may  be  computed  by? 


-a 


tan|(=  -F 


d'-y 


,  _  -sin  K  -cos  K  _  _JL 


■-c. 


2  2 
a  +b 
o  o 


^i 


’ll 


-3^2+^)+  ^(a^bg) 

2  ~2 
a  +  b 
o  o 

-aoW-b2)-bo(a2*bl) 


a2ab2 
o  o 


cot  (A+K)  = 


ao(bo  “l  -  ao  hi1  -  bo(ao  b2  -  bo 


-bo(bo  “l  -  ao  V  -  ao(ao  b2  -  bo  a2) 

a  (bQ  ^  -  a  b^)  -  b  (a  b  -  b  a  ) 
h' - *  cos  (Mtf -  <d*  - 

-bQ(bo  a1  -  aQ  b^  -  aQ(a0  b2  -  bQ  ag)  , 

- a  alb  (A +4 -  (4'-V> 


Or  introducing  and 


cot  (A+K) 


"Mi  -  b: 


a 

O 


+  a. 


»V  ”  al  “b0^i  “b2 


And 


(“%>\i-al)(d,-yc,). 
*  +sin(A+i() 

f-  (ao^i'*  ^2)(d,"yc,) 

*  +cos  (A+lO 


(-b0ni  -bi»d'  -  ?c2 

+CQS  (A+>0 

(-bo^l  -  V(d'  - 

- asii  a*K) - 


1 


(15) 
cont  * d 


(16) 
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i 


d'-y  A, 

C. 

.  'h' 


ye$y-p  ;|>$b2+  QVhti tt  A  -  Cc^*  y^^coi  '.  ‘!R~  — 
J "  ■  h* 


d*  »  (d  -  y^)  +  y^ 


d  =  d  sin  v 


sin  v  »  ~* 
h 


-  d  tan  | 


Ax  -  -  oos  £+  (y-Q  -  y^  sin  £ 

Ay  -  (y^-  y^)  cos<+  sin  M 

The  formulas  (lli)  again  express  the  relation  between  the  standard 
coordinates  £  ,  77  and  the  plate  coordinates  x,y  in  any  rectangular  plate 
coordinate  system.  The  unknowns  of  the  solution  are  the  so-called  plate 
constants  a^,  b_,  a,,,  a  ,  b  and  a_,  b_,  c2,  ao",  b  '.  Formulas  based  on 
plate  constants1  were  uied  in  earlier  worK  on  °ballistic  cameras  and  re¬ 
lated  problems  and  are  described  in  various  BRL  Reports.* 

*  1®  R.  S»  Zug,  "A  Photometerorohic  Method  for  Bomb  Ballistics  and  for 

Measurements  of  the  Flight  Fterformance  of  Aircraft,**  BRL  Report  279,  p.3£. 

2.  T.  E.  Sterne,  **The  Precise  Determination  of  the  Position  of  a  Point 

in  Space  from  Photographs  taken  at  Two  Ground  Stations 'J  BRL  Report  273,p.  111. 

3.  R.  S.  Zug,  "High  Altitude  Range  Bombinp  by  the  Aberdeen  Bombing  Mission 
Using  Ballistic  Cameras",  LSD-Report  1,  pp.  UO,  U3. 

li.  L.  G.  deBey  and  R.  S.  Zug,  "History  of  th<-  Ballistic  Camera  and  Range 
Bombing  Instrumentation",  BRL  Report  £77,  p.  2. 

£.  R.  d’E  Atkinson,  "A  Shortened  Method  of  Star-reductions  for  Calibrating 
Ballistic  Cameras,  BRL  Report  £93,  p.  12. 

6.  Virginia  Farquar,  "A  Short  Method  for  Reduction  of  Ballistic  Camera 
Observations",  BRL  Technical  Note  £0,  p.  2. 

7.  Dirk  Reuyl,  "Full-Scale  Free-Flight  Ballistic  Measurements  of  Aircraft 
Rockets  by  Optical  Methods",  BRL  Report  676,  p.  1*0. 


*U6) 
cc.;,‘t  *  d 


'X 

o  r  ■/ 


-sin  A  *  >  (o2+  dots  A 
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A  comparison  between  formulas  (lit) ,  (10a)  and  (10b) ,  which  obviously 
express  the  same  relations  shows  that  six  unknowns  are  present  in  formulas 
(10a)  and  (lOb)  and  ten  unknowns  in  formula  llu  We  have  already  seen  that 
six  unknowns  are  necessary  and  sufficient  to  solve  the  problem  geometri¬ 
cally. 


Consequently,  the  ten  unknowns  in  formual  (lU)  are  not  all  inde¬ 
pendent  parameters,  hit' four  additional  condition  equations  must  exist 
between  at  least  some  of  the  unknowns.*  A  first  examination  of  the 
terms  expressing  the  plate  constants  in  formulas  (l5)  shows  that 

1;  a  a  a 
o  o 

2)  b  5  b 
o"o 

The  other  two  condition  equations  can  hardly  be  eliminated  from  the 
given  relations.  However,  it  is  possible  to  write  the  equations  (10a) 
in  the  following  arrangement? 

y( -cos  <2.  cos  v  cos  A  +  sinA£  sin  A) +x( -cos  Ad  sin  A  -  sin  K cos  v  cos  A)  + 

4y(cos*lcos  v  cos  A-sin<sin  A)*^x(cos  Atsin  A+sinAdcos  v  cos  A)-d  sin  v  cos  A 
-y  cos  £  sin  v-x  sinAdsin  v+4y  cos  Ad  sin  v+  4  x  sinAdsin  v+d  cos  v 


y(— cos Ad cos  v  sin  A-sinA^cos  A)+x(+cosAtf  cos  A-sinJfcos  w  sin  A)  + 


4y(cosA£cos  v  sin  A+sin^cos  A)+4x(-cos  Adcos  A+sinAdcos  v  sin  A)-d  sin  v  sin  A 
-y  cos K.  sin  v-x  sinAdsin  v  +  Ay  cos K  sin  v+4x  sinAdsin  v  +  d  cos  v 


$  -i  *  *  bi  r  *  ci 

a  x  +  b  y  +  1 
o  o 


a;  *  ♦  b8  y  *  V 

a  x  +  b  y  +  1 
o  o  * 

where 

-  cosAdsin  A  sec  v  -  sinAdcos  A 
*1 - \ - 


*  See  ERL  Memorandum  Report  U76  by  S.  T.  Zaroodny,  "On  the  Use  of 
Least  Squares  in  the  Determination  of  Plate-  Constants." 


29 


h)  aQ2  +  a12  +  a22  -  bQ2  -  “  b2  =  0 

Using  formula  10^,  the  plate  constants  are  given  by  the  following 
expressions: 
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The  condition  equations  may  be  proved  by  substituting  the  expressions 
for  the  plate  constants  obtained  either  from  formulas  l£,  17a  or  17b. 

Formulas  lii  with  the  ten  plate  constants  and  the  four  condition  equa¬ 
tions  (formulas  19)  represent  now  a  system  of  six  independent  parameters. 
Now,  all  three  results ,  obtained  from  the  two  methods  are  consistent  and 
in  agreement  with  the  theory  of  the  orientation  problem. 

All  three  formula  systems  (10a,  ,10b and  lit  /  19)  express  in  rigorous 
geometrical  or  analytical  terms  the  relation  between  the  standard  coordi¬ 
nates  and  the  corresponding  plate  coordinates.  Each  of  the  three  sys¬ 
tems  allows  the  computation  of  the  plate  orientation  from  three  stars. 

In  general,  there  will  be  more  stars  used  than  are  necessary  for  a  unique 
solution.  In  such  a  case,  the  result  must  be  determined  bjr  a  least  squares 
solution,  which  may  be  based  on  any  of  the  three  systems  of  formulas . 

3.  The  rigorous  least  squares  adjustment 

The  formulas  10a,  10b,  and  1)4,  the  last  in  combination  with  the 
condition  equations  19,  represent  the  rigorous  mathematical  relation 
between  the  plate  coordinate^  measured  in  any  plane  rectangular  system 
on  the  plate  and  denoted  by  x  and  y,  and  the  standard  coordinates  £  and 
Vf  computed  as  plane  rectangular  coordinates  in  a  plane  tangent  at  a  cer¬ 
tain  point  on  the  unit  sphere  whose  center  is  at  the  point  of  projection. 


31 


All  three  systems  are  therefore  suitable  for  establishing  the  observa¬ 
tion  equations.  Jpfe  will  consider  these  possibilities? 

a)  Projection  method,  based  on  the  formulas  10a 

The  observations  of  the  plate  coordinates  with  reference 
to  a  fiducial  marks-system  are  denoted  by  i.  and  ,  corresponding  to  the 
x  and  y  axes.  The  observational  (residual)  errors  of  these  observations 
are  denoted  by  V  and  v'  j  hence,  x  =  A  +  V  and  y  *  T  Ax 

and  A  y  are  the  cuordinauea  of  the  principal  point  (P)  with  reference  to 
the  fiducial  mark  system  (x,y)  whose  origin  is  in  0. 

d  =  principal  distance 

A  «=  azimuth  of  the  optical  axis 

v  »  tilt  of  the  optical  axis 

fC  =  swing  angle  of  the  fiducial  marks  system 


From  formula  10a  we  obtain 


(y-dyXcps/gcps  v  cos  A°sin/t*sin)^-(x-43fr (sin/^cos  v  cos  A+cos/g slnA)+  d  sin  v  cos  A  ^p 
(y-Ay)  cos^sin  v  +  (x  tAx)  sin^sin  v  -  d  cos  v 

(20) 


(y-AyXcos^ cos  v  sin  A+sin/gcosfl)fr(x"dxXgiP^cos  v  sin  A-cos/f'cos  iQ-td  sin  v  sin  A^p, 
(y-dy)  cos/^sin  v  +  (x-iAx)  sintfsin  v  -  d  cos  v 


From  the. Taylor  expansion  for  the  right  side  of  these  equations  we  have, 
neglecting  terms  of  second  and  higher  order; 


dF 

3v 


4  u  ♦  fs4*- ♦  lkd  + 


3F 


3F 


dF 


dF 


=-AAX  +ZZ-AAV  +  —  V  +  — 

^  aEc  ay 
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♦?> 


3F»  . 

+  ~Av  + 
3  v 


3F*.  3F»  A 

H4  ^  +  a5kAAx  + 


3F 


y 

as? 


3F» 
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or  . 

V  +  -  V* 
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We  have 


4 c-So 


’Jc'l.  - Ar) 

where  4"  cJ  ^?c  are  standard  coordinates  computed  from  reference  data. 
Hence 

**  —  v-  “  *'  -  3ft  A A  +5T  4  V*  5*4*  +  3d  4  d+3^  y'4£ 


3F5  dF  ,  3F«A  »  SF»  3F«  Ai(/  aFf>  j  3F'  t  3F»  „  .  A  ~ 
—  v~  —  v  -  gc  A  I*s«  V*  -^,AK*  5(4  d*a£c*4  y- 


We  substitute  now 
A  =  Ag  +  4  A 

v  =  v  +  4  v 
o 

+  A< 

o 

d  =  d  +  4  d 
o 

A  X  -4xq  +  ^x‘,A7  =^y0  + 

5  ■  l  +  v 

y  *  Jt*  +  v® 

and  introduce  the  following  terms s 

A  (A*  -4y  )  *  B(A  -4xQ)  +  C*d0  s 

£  a  - - - - - - -  =  — 

E  (J>  -Ay0)  +  EU.-4x0)  +  F*dQ  u 
For  the  computation  of  approximation  values  see  page  50. 


(22) 


33 


A'U'-dyJ  +  B '(A  -^x J  +  C*  *d 


Vo  =  D(A'-^y  J  +  E  (i,  -4x  )  +  F*  d 


o 

a  — 

o  u 


where 


cos 

A  -sin<. 
o  o 

sin  A 

o 

A' 

=  cos  AC 

o 

COS  V 

0 

sin  A  + 
o 

sin  /t*  cos  A^ 
0  0 

cos 

A  +COS  /m 

sin  k 

B* 

=  sin^tf* 

COS  V 

sin  A  - 

cos/C  cos  A 

o  o 

0 

o- 

0 

o 

0  0 

o 

C« 

=  sin  v 

o 

sin  A 

o 

G  ■  sin  v-  cos  A 
o 

D  -  cos  K „  sin  v 
o  o 

E  sin  K  sin  v . 

o  o 

F  =  -  cos  v 


From  formulas  20  and  21  we  obtain  after  suitable  transformations 
the  observation  equations? 


■B+t  £  -A+  £  D 

»  O  IQ 


u 


V  + 


v-  =  ->?04A 

-[(t+4o}  cos  V^o  sin  Ao]^V 


-A+  t  D  -B+  ij  E  ~ 

-  “V —  -U'-4y0)  - j  At 


u 


-M  oB 


u 


'AA~x. 


A  Ay 


(23) 


°-$oF 


u 


A  d 


with  the  weight  p  ■ 


(-■MoE)2 


♦  '(-A^0E)2 


The  entirely  analytical  computation  of  some  of  the  coefficients;  asks  for  rather 
complex  transformations.  The  result,  however,  may  be  obtained  oirectly  by 
using  vector  analysis.  Such  a  solution  is  shown  in  BRL  Report  No,  785, 

R.  Schmid,  ^Error  Theory  of  Intersection  Photogrammetry.'* 
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-8  ’*nj  -k'*n  oE 
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V* 
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V'=*  +^q<Aa 


-[ 


(l+y?  )  sin  A  +n  cos  A 
•  o  0  10  o 


‘J 


-  (i-4x  )  — — 

1  o'  u 


-A«  +  ^D  -B 


-  U'-4y0> 


1  +n0*  1 


Ate 


-B,*n0E 


U 

=A»+7JqD 
u 


0  y 


(2T  cont'd) 


c'-n0F 


u 


A  d 


-AV 


with  the  weight  p8= 


u 


(-B'  +  rj  K)2«(-A  >+^  7^ 


The  number  of  unknowns  may  be  reduced  by  eliminating^  A.  We  divide  for 
this  purpose  each  of  the  observation  equations  by  the  factor  which  is  com¬ 
bined  with  A  A  thus  obtaining  reduced  observation  equations „ 

We  introduce  furthermore; 


+B-  tt  E 
*  0 

a*  = 

E 

t 

s 

♦M0d 
—  - 

KJ  B 

-A«+  n  d 

*  0 

w 

-Cf$oF 

f1  ?  B 

s 

+G'-r}? 

t 

s 

-  *£_ 

L»  = 

♦ 

no 
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Formula  23  may  be  written  as  follows? 


p=a*v>b.v  »JA+ 

‘  l+£  2 

A  v- 

(i.-4xo)b-U'-Ay  )a 

„  O  p  0  o 

L  n° 

J 

+a*44x+b»44y+c»4  d-L 

p’sa’v-b'  v*=AA- 

i:V 

4  v- 

(A-4xo)b'-(A'-4yo)a 

fe  0  1°  o 

o  0 

- 

+a'  •  44x+b*  •  44  y+c'  • 

4  d-L' 

Af£ 


(230 


or: 

p=  4A+a4V+p4£  +a»AA  x+b°  44  y+c»4d-L 
p'=4  A+a'4  v+P'4  Ai  +a'  •  44  x+b'44  y+cU  d-L* 


and  with  the  help  of  the  first  normal  equations ,  we  obtain  the  reduced 
observation  equations: 


p-  (A)4  v  +  (B)4*L+  (C)44X  +  (D)44y  +  (E)Ad  -  (L) 
pi-  (a04  v  +  (B04K+  (C» )  Aix  +  (D')44y  +  (E»)4  i  -  (LO 


(2U) 


(A) 

=  (a 

[a+a*J\  ' 

"  ~2rT  )> 

(AO 

m(a* 

[a+a'J  \ 
2n  J 

(B) 

JS1S3 ) . 

2n  /» 

(bO 

-(*' 

JM9) 

2n  >/ 

(c) 

s  (a 

[a+ag\ 

2n  /' 

(CO 

[a+aQv 

J 

(D) 

“  ^b 

^ 

(DO 

=^b* 

■  Cb+b'jA 
5n"  / 

(E) 

=  (c 

[c+c'3\  . 

“  “?ir  y  > 

(BO 

=(c 

(L> 

-(* 

cl+l  a  \ . 

’  “sr  /' 

(lO 

4' 
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and 


A  A-  -[2^4  v 


Fu  *  g 


(2^) 


v  =  V-p-^tei. 

ab’-a'b 


(26) 


v» 


ap'-a'p 

ab«-a»b 


The  reduced  observation  equations  (2i|)  have  different  weights  due  to  the 
fact  that  in  each  of  the  equations  there  is  more  than  one  observed  quantity. 
The  weights  ares 


P 


>»  * 


+  b5< 


(27) 


The  normal  equations  ares 

|p(AA)]4v  +  [p(AB)|4£  +  [p(AC)]  44  x  +[p(AD)]^4y  +  [p(AE)J/Jd-  |p(ALj] ■  0 
[p(BB)|^/e:  +  [p(BC)]4*i  +[p(BD)]j4  7  +[p(BE)]4i-  [p(BL)]  =  0 
[p(CC)]  44  x  +  [p(CD)j44  y  +|^(CE)]4d-  |p(CL)]=  0 
[p(DD)]44y  +|^(r)E)]4dl-  [p(DL)]  =  0 

Jp(ee)]4A-  [p(el)]-  o 
[p(LL)]»  0 


m 


where 

[p<U)|  stands  for  |p(A)(A)  +  p,(A,)(A,)j  and  the  other  terms  in  parentheses 
are  denoted  correspondingly. 
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After  A  v,  .4/fcf,  A  A  x,  /Jy  and  are  determined  ,  is  obtained  from 
equation  25.  The  residuals  p  and  p*  are  computed  from  the  reduced  obser¬ 
vation  equations  2b  and  the  residuals  of  the  original  plate  measurements 
are  obtained  from  formulas  26, 

The  final  result  is  obtained  by  applying  the  formulas  (23). 

Checks  for  the  computations  ares 

|  [p]  +  [p']"  0  [PPP  +  P’P'P1]]  “  [vv]  +  -[W3  -jp«ll»5]  (tf) 

l 

The  main  check  is  made  by  means  of  formulas  10a.  If  the  final 
elements  of  orientation  and  the  corrected  observations  x  and  y  are  intro¬ 
duced,  the  computed  f  and  ^values  must  be  in  complete  agreement  with 
the  originally  obtained  standard  coordinates  $  and  7J  respectively. 

In  case  of  residuals  which  are  caused  by  the  neglected0  second  order 
terms,  the  solution  must  be  repeated  by  using  the  obtained  results  as 
new  approximations.  It  will  be  sufficient  in  such  a  case  to  use  the 
former  coefficients  in  the  normal  equation  system  and  to  limit  the  addi¬ 
tional  computations  to  the  computing  of  new  L  and  L*  terms. 

The  mean  error  of  a  measured  plate  coordinate  may  be  determined  by! 


.(30) 


and  the  mean  errors  of  the  unknowns  are! 


m  =  + 
v  — 


m 


m 


\/aaT4" 


I  m  ■  +  - 

AX  -vrccT 


m 


“d  - 


\/  EE*U 


(31) 


V3  ! 


m 


m  a  + 


m 


Vbb-u  *  vto*u 


For  the  reductions  of  the  normal  equations  the  M  modernized”  Gaussian 
Algorithm  is  suggested.  The  method  is  closely  related  to  the  Doolittle 
method  and  to  Reicheneder!s  method.  However,  it  is  believed  that  the 
modernized  Gaussian  Algorithm  combines  a  minimum  of  recording  work  with 
a  maximum  oi  protection  against  computing  errors.  Sees  Der^njodemisd erte 
Gauss !sche  Algorithmus  zur  Auflosung  von  Norma Igleichungen,  by  H.  Wolf  in 
Zeitschrift  f£Lr  Vermes  sungswesen  ^11  November  1950. 
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m.  may  be  determined  from  formula  (25)  by  combining  the  individual  mean 
errors.  However,  the  mean  errors  are  not  independent  and  therefore  can¬ 
not  be  propagated  according  to  the  Gaussian  law  of  errors.  The  values 
AA°4 >  etc.,  have  here  the  conventional  meaning  of  the  Gaussian  reduction 
method.® 

b9  Projection  method,  based  on  the  formula  10b. 

The  observations  of  the  plate  coordinates  with  reference  to  a 
fiducial  marks  system  are  again  denoted  by  1  and  1*  corresponding  to  the 
x  and  y  axes_.  The  observational  errors  of  these  observations  are  V  and 
y’  Hence,  x  =  A  +  v  and  y  =  JLl  +  v* .  A x  and^y  are  the  coordinates_of 
the  principal  point  (P)  with  reference  to  the  fiducial  marks  system  (x,y) 
whose  origin  is  in  0.  d  denotes  the  principal  distance,  a  =  £  -tilt  angle 
and  co  ~T]  -tilt  angle  respectively.  fC  denotes  the  swing  angle  of  the  fiducial 
marks  system.  From  (10b)  we  obtains 


(y°Ay)(sin£cos  a-cos<sin  oo  sin  a)+(x-4x)  (cos  <cos  a+sin£sin  co  sin  a)+d  cos  cosin 
(^tiy)  (-cos£  sin  as  cos  a-sinAS  sin  a)+(2-4x)  (sin^Csin  co  cos  a-cosj(£sin  a)+dcoscocos  a 


(32) 


(y-4y)  cos^cos  co  -  (x-^x)  sin£cos  co  +  sin  co 


■G* 


(y-Ay)  (“Cos/Usin  co  cos  a-sin£sin  a)+(x-Ax)  (sin^sfn  co  cos  a-cosjK  sin  a)+d  cos  co  cos  a 


From  the  Taylor  expansion  for  the  right  side  of  the  above  equations, 
we  have,  neglecting  terms  of  second  and  higher  orders 

*  -  *  §  * «  Sx—  1  v  *  | v' 


V -%  *S"°*  S’" -  *  S’ " d*  E "" -  l'v*|’  v 


and  with  ^  c  -50 

r/  c  -rja  •  Ar) 

where  and^c  are  standard  coordinates  computed  from  reference  data., 


(33) 


Refs  A  Text  Book  on  the 
N.Y.,  1911;  ors  Handbuch 
Jordan  Eggert,  volume  I, 


Method  of  Least  Squares  by  Mansfield  Merriman, 
dar  Vermessungs^JCunde-Awsgleichsrechnung  by 
Stuttgart,  19U8. 
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.a  v. M 

Edc  3y 


v>=  |Ua+ 


3G 

dco 


Zco+ 


8**  i  « 


.  3G 

dt  " 


3G 

diy 


AA 


y-  A  ^ 


Ml  v-  Mf 

3k  3y 


v’=  a+  a* 


da 


MV:+  S’*d+  301 


dK 


Sd 


ddx 


§y  A*y- 


We  now  substitute 


a  ■  aQ  +  4  a 

co  =  co  +  A  co 
o 

iC  +  4AZ- 

o 

d  =  d  +  A  d 
o 

4 x  ■  4x  + A A  x 

o 

4y  »'^y  +  AA  y 
x  -  1  +  v 
y  -  2'  +  v" 

and  introduce  the  following  terms! 


(3U) 


A(A‘  +  B(jt  -4xq)  +  Odo  a 

i  a - - - — -  ■  — 

70  DU'  -JyQ)  +  EU  -4xo)  +  ?.do  u 

A'  (A1  -4y  )  +  B'(i  -Zx0)  +  C*  *dQ 

T9  *  -  °  — 

1  0  D(l'  -4y0)  +  E (1  -iixQ)  +  F*do  u 


(35) 


where! 

A  =  sin£0  cos  aQ-  cosjtfQ  sin  ooQ  sin  aQ 

B  “  coS/K  cos  4  +  siniK„  sin  00  sin  a 
o  o  o  o  o 

C  ®  cos  00  sin  a 
o  o 


A'  =  cos  cos  a)Q 

B*  «-sinA£  cos  co 
o  0 

C*  =>  sin  a) 


Uo 


D  =  -  cos  K  sin  at  cos  a  =  sin  AS  sin  a 
o  o  o  o  o 


E  *  sin  £  sin  as  cos  a  -  co a  £  sin  a 
o  o  o  o  o 


F  =  cos  a>  cos  a 
o  o 


From  the  formula  (32)  we  obtain  the  observation  equations  after  suitable 
transformations* 


B+  $  E  -A+i  D 
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+d  ,-£02)'4a 
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o'  u 


-A+£  D  ~B+t  e" 
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-A.  f0D 
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4  4  7 
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■4^  with  the  weight  p  = 
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(-B*^oE)'.(-A*f0D )' 
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-B*+^oE  “A'+  ^ 
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V*— S-5-  V'  m*rloto‘ia 
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See  remark  on  page  3U 
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A  A  X 


-B'+>?oE 

U 

u 

C'->?oF 
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AA  y 


A  d 


-+n  with  the  weight  p' 


u 


(-B>  +  /7oE)2+C-A‘+^bD)2 


The  number  of  unknowns  may  be  reduced  by  elim±nating4o.  For  this  purpose 
we  divide  each  of  the  observation  equations  by  the  factor  which  is  com¬ 
bined  with  A  a  thus  obtaining  reduced  observation  equations. 

We  introduce  furthermore: 
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-B  +  F  E 
i  o 

3  ■■■■ 

u  +fe  *s 

•  i  o 


-A 

U  *$o-3 
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+C  4  F 
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Formula  36  may  now  be  written  as  follows; 

T)  £  cos  a  -  n  sin  a 
10*0  0  to  o 


p=av+bv*=^<i+ 


Cl 


Aon*- 


Ci-4  XQ)b-U:*  - 


4y0)a 


AK+ 


(36») 


a  -44x+b  a 4  y+cA&  -  L 


4jtr+ 


p  •  «s  a  *  y  *b  •  v  *  ■  ^  a+ 


bo?, 


sin  a 


cos  a  + 
o 


o  _l 


Av* 


Uftd^b’-Ci'  -*4  y  )a 


a8  <A4x+b'  44y+cMd '*»  L1 


where  L 


Ho 


L>  = 

no$c 


The  weight  of  the  observation  equations  is  p  *  -5 - 

aT  +  b^ 


p»  = 


1 

T~ — ? 
a'^+  bf  * 


(37) 


The  solution  now  follows  the  procedure  outlined  in  the  preceding  paragraph. 

c)  Analytical  Method,  based  on  formulas  lit  and  19 
Substituting  aQ  a  a0*  and  bQ  a.  bQ'#  formulas  lh  become 


and 


V*  V  *cl-  »o?*-  »0?7-<  -0 


a2x  +  b2y  +  c2  -  aQ?7  x  -  bQ7jy  -  7J  =  0 


(38) 


We  again  introduce  x  »  1  +  v  and  y  =  1*  +  y’  and  approximation 
values  of  the  plate  constants  as  follows 


o  . 

a1  = 


a2  «  a2°  *  A  a2 


bl  "  bl°  +4bi  * 


b2  “  b2°  +  4 b2 


U3 


(39) 


°1  «  °1°  ♦4c] 


C2  ®  C2  +  A  c2 


a  *  a  +  Aa.  . 

0  0  O  9 


b  »  b  °  +  4b  , 

o  o  o  * 


W.e  obtain  from  formulas  38,  by  the  Taylor  series,  neglecting  terms  of 
second  order  and  highers 

p=X  A a1+  A  b^+  A C;L-i^  A  f  +la1°+t,b-L0+c1°”Jl^  ao°-l'^  b0°  " 

ao°-a1°)4V'  ($  bo°-b1°) 


p*«X4  ag+Jl’  4  b2+^  c2-X?7^  aQ-X'  77  ^b^^^^'b^Cj^E^a^-i^  b 


W) 


O  ,0> 


pai-ila^+t'd  b^  A  c^-1^  c-dao-X1  ^  c4bQ-L,  with  weight  p 
p»=i4  a2+£‘4  b2+4  c2-l?7c4  ao-l*7)  c4bd“L’ ,  with  weight  p' 


p’(fao“.aio)^(fboo-b1o)2  5  r  <  v0°V>z*(,>  VV>Z 

The  introduction  of  weighting  factors  becomes  necessary  because  in  each 
of  the  above  derived  observation  equations  more  than  one  observation  and 
correspondingly  more  than  one  residual  appear. 

-L  -  1  *1*  £V  *  cl°  -  lt  cao°'£l  f  cbo°  -  t  c 
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“L*=  la°  +  !'b,°+  c°-  in  a  °-  Vn  b  °  -n 
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where  the  computed  standard  coordinates  are  denoted  by  f  c  and^l  c’ 


uu 


The  condition  equations  (19)  must  be  satisfied  by  the  unknowns.  The  first 
two,  namely  a  a  a^5  and  bQ  a  bQ?  have  already  been  introduced  into  the  equa¬ 
tions  (38).  The  ° third  and  fourth  conditions  ares 


3.  a  b  +  a.b.  +  a0b0  =  0 
o  o  11  d  d 

U.  a/  *  *  a/  -  bo2  -  bx2  -  b/  -  0 

By  means  of  the  substitutions  made  in  formula  38  and  the  Taylor  series, 
neglecting  terms  of  second  and  higher  order,  we  obtains 

bl0<ial  +  al°4bl  *  b2°4a2  *  a2°4 b2  *  bo°4  *o*V4bo  *  h  ‘  0 


'*ere  h  ’  aobo  *  hh  *  “2b2 


and 

O  i  O  >  O  a 

ao  4  V  *1  4*1  *  a2  4a2 


bo°4  bo  -  bl°4bl  *  b2°4  b2  *  \  - 


where  X, 


2  2  2  2  2  2  * 
ao  +al  +a2  “b0  “bl  ”b2 


After  suitable  transformation,  we  have 
lAa^+l'ATxf  *A  a2(ia»+g«b')+4  b2(£«a*-£b«  )+4  a^c'+j^d* )+  ^ bQ (Jt » c * -jtd *  )♦ 

£e'  +  £'f»  (U2) 


If  the  approximation  values  a  °,  a_°,  a_  ,  b  ,  bn  and  bp  are  chosen 
In  such  a  way  that  they  satisfy  the  condition  equations  3^ and  U,  both 
X.,  and  X-  become  zero.  At  the  end  of  this  paragraph  a  "procedure  for 
obtaining  such  approximation  values  is  shown.  In  routine  reductions, 
where  usually  the  preceding  result  will  be  introduced  as  the  approxima¬ 
tion,  the  condition  equations  will  in  general  be  sufficiently  satisfied. 


where 


o 


d» 


e* 


f » 


-  0  \ 

Ki  ”  n.  k2 


+b-,  °X„  -  a^° 


1  2 


w 

*1 


2  2 

O  .  .0 

+  V 


m 


For  »  0  and  =  0,  e'  =  0  and  f 1  =  0  and  we  obtain 

4  a,  =  a'4  a0  -  b'4  b_  +  c*4  b 

x  2  2  o  o 

m 

■4  b,  “  b'^a.  +  a'4b„  +  d’ 4  a  +  c'4b 

1  <?  2  O  O 

Substituting  h2  with  e’  =  0  and  f1  =  0  in  the  observation  equations  1*0, 

we  obtain  the  observation  equations*  1 

Ac^la'+l'b’)^ a2+a>a'-lb0^b2+(lc«+Jttd'-Jt^c)4ao+(l'c»-id>-l'^  ^b^t-p 

(U5) 

v  4b2-  £7c  4-v  *-xrlc  ^VL,ap‘ 


1*6 


4  c2  +  l  4  a2+ 


or 


+  a*da2  +  p^b2  +  Y A*0+S A  bQ-  L  =  p 
A  c2  +  fc^a2+  VA\-lT)  c^aQ-  L*  -  p» 


(U6) 


where 


a  =  flU'  +  jt’b') 

P  -  (i'a'  -  lb') 

T*  j£(c'-£c)  +  ?  <f-[£'(c>  ~  §c)  -  £d|] 

Before  we  form  normal  equations  we  eliminate  the  unknowns  A  c,  and  A  C£, 
Thus  we  obtain  the  reduced  observation  equations!  & 
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{&--)4a0+  (£•-  -  )4b0-(tn  - - )  a  -  (i ,*/l  -  - -  )Ab  -(L» - 
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) 


or 

A  Aa^-ft  B^b2  ■*■  C^aQ  +  DAbQ  -  (L)  =  p  with  the  weight  p 

<U8) 

A'Aa2+  B*  A  b2  +  C?A  aQ  +  D'A  bQ  -  (L1)  =  p”  with  the  weight  p* 

The  corresponding  normal  equations  ares 

JpAAjd  a2  +  (pABjA  b2  +  [pACj ^  aQ  +  JpAD] A  bQ  -  £pA(L)J  =  0 

[PBB>  b2  +  [pBc]  A  aQ  +  [pB^Z  bo  -  [pB(L)]  «  0 

[pGC]4  ao  +  [pCDj  A  b0  -  QpC(L)]  =  0 

[pD§  AbQ  -  £pD(L)J  =»  0 

[p(L)(L)|-  0 

For  these  equations  jpAAj  stands  for  [pM  +  p'A!a]  and  the  other  terms 
in  parentheses  are  denoted  correspondingly. 
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After  the  unknowns  A  a^,t  Ab^tAa  and  A  bQ  are  determined*' we  compute  A 
and  A  Cj  by  solving  uhe  equations:  0 


r«]  [fl  cn  m  [13 

-4a*  -  -ib»--Za  -  —  4b  +  — 
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(50) 


and  from  the  equations  we  compute  A  and  A  b^ . 

The  final  plate  constants  are  obtained  from  the  equations  (39) . 

As  a  first  check  these  values  are  introduced  in  the  condition  equations 
(formula  19),  which  must  be  satisfied.  In  case  the  condition  equations 
are  not  sufficiently  satisfied  the  solution  must  be  repeated  by  using  the 
results  as  new  approximations.  (  Compare  with  p.  38  )  In  such  a  case 
it  is  advisable  to  compute  the  su  and  b^  values  with  formulas  57  in  order  to 
obtain  approximation  values  which  satisfy  the  condition  equations.  The  resi 
duals  p  and  p*  are  computed  with  formulas  i+8  and  the  final  residuals:  of  the 
plate  measurements  v  and  v'  are  obtained  from  (1+0)  as  follows: 

P(^bo0-b2°)  -p'ttb/-^0) _ 

ao°“al°)  bo°“b2°^"^  ao°“a2°^  bo°-bl°) 


(?D 


v 


f"<  Vo°-al°)  -  0^ao°V) 


ao°-=l°>  1 1  bo°V> 


A  check  is  obtained  by 

D3  =  &>[]  “  0 

(  *  • 

and 

Cvv3  ♦  Cv'v'J  =*  CpppH  +  Crflp'pn  -  p(ll)  •  3  (52) 

The  finalcheck  is  made  by  means  of  formula  li+.  After  the  x  and  y  values 

are  pbtained  with  x  =  £<■+  v  and  y  ~  t1  +  V  the  corresponding  standard 

coordinates  4  and  ?7  must  be  in  complete  agreement  with  the  originally  com¬ 
puted  ^  and  7)  .  In  addition  the  unknowns  must  satisfy  the  condition 

equationS  19 •  C 


Compare  footnote  on  page  38. 
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The  mean  error  of  an  observed  plate  coordinate  Iss 


fvvQ  +  Cv'y'3 


m  = 


2n  -  6 

and  the  mean  errors  of  the  unknowns  may  be  computed  from? 


(S3) 


m 


jna  =  +  . . — 

2  \TaTJ 
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*  “  Vbb*3 
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(Si) 


ma  =  +  == 

0  ~4cc- 
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mb  *  +  . 

0  “VdD‘3 


The  values  AA*3>  BB®3*  etc.,  have  the  conventional  meaning  of  the 
Gaussian  reduction  method.  See  footnote  on  page  39. 

The  mean  errors  of  the  other  four  plate  constants  denoted  by  mal, 
m.  , ,  m  ,,  and  mc2  may  Be  obtained  from  the  formulas  Ui  and  £0»  The 
correct  computation  is  rather  complicated.  Since  the  individual  mean 
errors  of  the  unknowns  are  not  independent,  the  computation  must  follow 
such  lines  as  given  in  the  computation  of  the  mean  error  pf  a  function 
of  the  unknowns.*  For  this  purpose,  the  weight  of  these  functions  must 
be  determined.  If  the  corresponding  weights  are  denoted  by  Pal'  Pbl> 
P.t ,  and  P_9,  the  mean  errors  are? 


'clJ 


‘c2J 


mal 
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t  m\jT~ 

Vrbl 
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c2  - 


(55) 


Handbuch  der  Vermessungskunde  Ausgleichsrechnung  by  Jordan-Eggert, 
Volume  I  -  Stuttgart  19iiB  -  p.  99 
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All  three  least  squares  solutions  described  assume  that  the  reference 
data  -  i»e»,  the  computed  standard  coordinates,  are  free  of  errors.  The 
standard  coordinates  $  and  7J  are,  in  our  case,  functions  of  the  right  ascen¬ 
sion  and  declination  values  which,  in  turn, are  results  obtained  from  ob¬ 
servations  and  therefore  affected  by  uncertainties  which  may  be  expressed 
in  the  form  of  mean  errors.  If  the  propagated  mean  errors  of  the  standard 
coordinates  are  not  negligible  in  relation  to  the  expected  mean  errors  of 
the  plate  measurements,  it  becomes  necessary  to  make  allowance  for  correc¬ 
tions  to  the  computed  standard  coordinates.  In  such  a  case,  we  introduce 
$  m  %  c  +  v$  and  V  =7?c  +  *  respectively.  Such  a  step  calls  for  addi¬ 

tional  weighting  factors.  Unit  weight  of  an  individual  plate  measurement 
may  be  assumed  for  both  the  x  and  .y  measurements,  or  p  ■  p  •  k-  ■  1, 

—  —  x  ■  y  m2 

where  m  is  the  expected  mean  error  of  the  x  and  y  measurements.  If  the 
propagated  mean  errors  of  the  standard  coordinates  are ’denoted  by  mi  and 
the  corresponding  weights  will  bet  ^ 


vn 


k 

{mrj  *d)^ 


(56) 


The  fact  that  the  orientation  of  the  ^  system  does  not  necessarily 
correspond  to  the  orientation  of  the  x,y  system  must  be  considered  in 
computing  the  corrections  v  £  and  v^  , 

The  Computation  of  Approximation  Values  for  the  Unknowns . 


Approximation  values  for  the  .unknowns  are  needed  in  all  three  methods 
of  adjustment.  In  the  projection  method,  approximation  values  of  the 
orientation  elements,  and  in  the  analytical  method  approximation  values 
of  the  plate  constants  are  necessary.  The  latter,  in  addition,  should 
satisfy  the  existing  condition  equations  (19).  The  relation  between  the 
orientation  elements  and  the  plate  constants  and  their  inverse  transforma¬ 
tion  are  given  by  formulas  l£,  16  ,  I7a  and  17b, 

In  routine  reductions,  approximation  values  will  generally  be  avail¬ 
able  from  the  results  of  the  preceding  reduction.  In  the  case  that  none 
are  available,  the  following  procedure  may  be  followed?  With  five  stars 
well  distributed  over  the  plate,  the  plate  constants  are  determined  from 
formula  (38).  From  the  first  equation,  we  obtain  a^,  b^,  c^,  ao’  and  bQ' 

and  from  the  second  equation,  we  compute  b^,  Cgj  aQw  and  bow .  The  aQ’ 

and  aQH  values  and  the  b  1  and  b  w  values,  respectively,  are  combined,  by 
forming  the  arithmetic  averages? 
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a  1  +  a  11 
o  o 


a 


o 
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and 
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o 


»  +  b  " 
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It  is  obvious  that  these  approximations  will  not  satisfy  the  condition 
equations  (19).  As  a  next  step,  therefore,  we  compute  a  set  of  orienta¬ 
tion  elements  using  formulas  (l6).  These  values  may  betaken  directly 
as  approximation  values  for  the  adjustment  of  the  projection  method. 

If  the  measuring  camera  is  well  adjusted,  it  normally  will  be  sufficient 
to  take  zero  as  the  approximation  values  for  the  coordinates  of  the 
principal  point.  Approximation  values  for  the  plate  constants  may  now 
be  computed  by  applying  formulas  17,a  or  17b.  The  final  approximation 
plate  constants  should  be  checked  by  the  condition  equations  19,  before 
they  are  introduced  in  the  least  squares  adjustment  of  the  analytical 
method.  For  the  analytical  method  a  set  of  approximation  values,  con¬ 
sistent  with  the  condition  equations  19  may  be  directly  computed  with 
the  help  of  the  equations 


"a2b2 


-a  b 
o  o 


The  method  of  least  squares  adjustment  which  should  be  applied  depends  on  the 
particular  situation.  In  one  case  where  all  six  elements  of  orientation, 
are  unknown,  the  analytic  method  seems  to  be  advantageous,  because  the 
least  squares  adjustment  can  be  arranged  in  suoh  a  way  that  only  four 
normal  equations  need  be  reduced.  However,  there  may  be  the  case  in  which 
some  orientation  elements  are  known,  e„g„,  from  a  preceding  calibration 
procedures  thus  calling  only  for  the  determination  of  the  remaining  un¬ 
known  elements .  The  formulas  given  in  the  ieast  squares  methods  for  the 
projection  method,  are  then  preferred,  because  some  of  the  unknowns  will 
be  equal  to  zero.  If  phototheodolites  are  used,  where  the  approximation 
values  are  obtained  from  dial  readings,  the  projection  method  offers  a 
decided  advantage . 
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IV.  THE  MATHEMATICAL  ANALYSIS  OF  THE  SPATIAL  TRIANGULATION  PROBLEM 


The  theoretical  background  of  the  triangulation  problem  is  explained 
in  Chapter  II.  There  it  is  shown  that  we  have  to  deal  with  the  two  fol¬ 
lowing  triangulation  methods? 

l)  The  rigorous  solution  in  triangulating  a  point  from  independent 
fixed  photogramme trie  cameras  calls  for  resection  in  space. 

2}  After  ^calibration  of  the  photogrammetric  base  lineH  the  tri- 
angulation  may  be  obtained  by  spatial  intersection. 

The  latter  method  has  been  described  in  many  reports.  A  method  of 
least  squares  adjustment,  suitable  for  combining  two  or  more  photogram¬ 
metric  measuring  results  for  the  purpose  of  triangulation  by  intersection, 
is  outlined  in  BRL  Report  752*.  However,  it  must  be  understood  that  al¬ 
though  the  described  method  provides  a  rigorous  least  squares  solution  for 
triangulation  problems  based  on  theodolite  data,  in  our  case  it  is  only 
an  approximation  method  due  to  the  fact  that  the  angular  corrections  and 
not  the  corrections  to  the  original  plate  coordinate  measurements  are 
minimized . 

We  will,  in  this  report,  deal  only  with  the  resection  problem. 

From  the  spatial  coordinates  of  the  measuring  stations,  denoted  by  x£ , 
y  ,  Zj_,  which  may  be  available  in  any  coordinate  system,  the  slant 
distance  between  any  two  of  such  stations  can  be  computed.  Assuming  the 
three  camera  stations  A,  B,  and  C,  we  have  the  .slant  distances  a,  b, 
and  c  between  these  stations.  (Fig.  10)  From  the  least  squares  adjust-*- 
ment  of  the  plate  orientation,  we  have  for  each  station  a  set  of  unknowns, 

either  plate  constants  or  orienta¬ 
tion  elements,  depending  on  which 
method  was  applied.  Introducing 
the  measured  plate  coordinates  for 
an  additional  recorded  target  point, 
after  suitable  corrections  for  non- 
perpendicularity  of  the  comparator 
axes  and  corrections  for  lens  dis¬ 
tortion  are  applied,  the  standard 
coordinates  ^  and  for  such  a  tar¬ 
get  point  are  determined  with  the 
formulas  10a,  10b,  or  lh  and  the 
X  corresponding  direction  angles  may 
be  found  by  solving  the  equations 
1  or  2. 


BRL  Report  752  by  H.  Schmid.  Titles  Spatial  Triangulation  by  Least 
Squares  Adjustment  of  Conditioned  Observations. 
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From  equation  1  for  the  resection  methods 


cot  6 

tan  RA 

and  from  equation  2  for  the  intersection  methods 
tan  z  2  +£  2 

^  m 

tan  A  =  — 

-1? 

The  spatial  angles  a  (Pig.  li  and  10)  at  the  apex  of  the  pyramid  ABCP 
may  be  computed  by  the  following  formulas s 


cos  a 


i in  6  sin  &  +  cos  6  cos  cos  A  RA  . 

*r  T7  r  r  aD! 


cos  <J  =  sin  6  sin  6  +  cos  6  cos  6  cos  A  RA 


cos  a. 


sin  6.  sin  6  +  cos  6,  cos  6, ,  cos  A  RA, 

r  °r  °r  r  D^i 


The  subscript  Nr*1  indicates  that  the  angler  are  corrected  for  refraction. 

(See  the  remark  about  refraction  on  page  1$ .  However,  either  astronomical 
or  terrestrial  refraction  must  be  considered  depending  on  the  spatial  po¬ 
sition  of  the  target  point.)  A  difference  in  time  (bonverted  to  sidereal 
time)  of  the  orientation  exposures  at  the  different  stations  goes,  into 
the  solution  as  a  corresponding  correction  to  the  right  ascension  values. 

Our  problem  is  now  to  determine  the  coordinates  of  a  point  P(xp,  yp,  zp) 

(Fig.  10) ,  First  it  becomes  necessary  to  compute  the  length  of  the  sides 
of  the  pyramid  denoted  tiy  s  ,  s,  ,  s  ,  and  finally  the  coordinates  of  the 
point  by  the  intersection  aof  Dthr8e  spheres.  The  second  part  is  a  familiar 
problem  for  a  triangulation  procedure  using  length  measuring  methods,  such 
as  Doppler,  etc.  A  possible  solution  is  outlined  in  BRL  Report  7li8.»  Hence 
we  are  left  with  the  problem  to  determine  the  lengths  of  the  sides  of  a 
pyramid,  given  the  length  of  the  sides  of  the  base  triangle  and  the  three 
angles  at  the  apex  of  the  pyramid.  This  is  the  basic  problem  of  the  re¬ 
section  in  space.  For  the  three  triangular  faces  of  the  pyramid,  we  have 
the  equations s 

*  BRL  Report  7li8,  H.  Schmid.  Titles  A  Rectangular  Ellipsoidal  Coordinate 
System  for  Trajectory  Measurements.  pp«llv-17,  pp. 23-25. 
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2  2  2 

a  =  s,  +  s  “2  s,  s  cos  a, 
be  be  be 


2  2  2 

b  =■  s  +  s  -2  ss  cos  a 

a  c  a  c  ac 


(61) 


°2  =  Sa2  +  3b2  "  2  sa%  cos  aab 


The  further  treatment  of  these  equations  results  in  an  equation  of  the 
fourth  degree  which  is  not  especially  suitable  for  numerical  computations. 
In  view  of  the  fact  that  an  economical  computation  of  the  final  coordi¬ 
nates  by  the  intersection  of  three  spheres  calls  for  the  introduction 
of  approximation  values  of  the  spatial  position  of  P,  it  seems  advisable 
to  follow  the  following  procedures 


by  x 


V 


Suitable  approximation  values  for  the  target  position  denoted 
and  z  must  be  obtained.  Such  values  may  be  computed  by- 


using  the  intersection  method  for  any  two  stations.  -The  necessary  azimuth 
and  elevation  angles:  may  be  computed  from  the  declination  and  right  ascen¬ 
sion  values  with  the  sidereal  time  and  the  geographic  coordinates  of  the 
station  by  well-known  formulas  of  spherical  astronomy.# 


sin£  ■  sin  fi  sin  6  +  cos  fi  cos  6  cos  t 
cos  6  sin  t 


(62) 


sin  A  - 


cos  £ 


t  =  hour  angle  =  sidereal  time  ($)  -  Right  Ascension  (RA) 


With  the  approximation  values  for. the  coordinates  of  point  P  and  the 
coordinates  of  the  stations  the  corresponding  approximation  values  for 
the  sides  of  the  pyramids  (s^  )  may  be  computed. 

Applying  the  Taylor  series  to  the  formulas  (6l)  and  neglecting  all 
terms  of  second  and  higher  order,  we  obtain: 


s 

a 


+  A  3 


a 


ab  "  sb°  +  ^  sb  (63) 

s  =s  +  A  s 
c  c  c 

See  references  on  page  IS  and  formulas  at  the  end  of  this  chapter. 
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ac  ac  a  ca  ac  c  0  a  c  ac 


(s1_°-s  °cos  a,  )4s,  +  (s  °”S,°cos  a,  )^s 
d  c  be  b  c  b  be'  o 


o  o. 
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— - -  s,  s  cos  a. 

a  DC  uC 


or 


*1  *  *  *1  *  =b  '  °1 

*2  A sa  *  B2  A  sc  '  °2 

A,  4s.  +  B,  4s  s  C, 
J  b  j  c  3 


and  therefore 

Where 


AS  = 

B-a 

DAS. 

fit 

a 

|D| 

Asb“ 

IXI 

|DI 

AS  ■= 

|D-ct 

DAabl 

lb.l 

VJ 

1b| 


CiBiO 

e2  o 


C3  A3  B3 


Ai  ci° 
A2  °2  B2 
0  C3B3 


Al0!0! 

A2  0  C2 
0  A3  C3 

AlBl  ° 
A2  0  b2 

0  a3  b3 


(6U) 


The  final  length  of  the  sides  of  the  pyramid  are  determined  by  the 
formulas  (63) »  The  determination . of  the  final  coordinates  now  follows 
the  procedure  outlined  in  BRL  Report  7U8.* 

During  the  orientation  and  the  triangulation  computations,  there  is 
sometimes  a  need  for  direct  or  Inverse  transformation  of  the  elements  of 
orientation  between  the  Right-Ascension-Declination  system  and  the  Azimuth' 
Elevation  System.  ■ 

The  following  formulas  may  be  used  for  this  transformation; 
d  »  principal  distance 

<=  Coordinates  of  the  principal  point 

=  Swing  angle  of  the  plate  coordinate  system 

v  =■  zenith-distance  of  the  plate  perpendicular 

A  *  Astronomical  azimuth  of  the  plate  perpendicular 
counted  clockwise  from  the  south 

RA  =  Right  Ascension  of  the  plate  perpendicular 
counted  counter-clockwise 

Declination  of  the  plate  perpendicular 

Standard  coordinates  of  a  star  in  the  azimuth- 
elevation  system 

Standard  coordinates  of  a  star  in  the  RA-declination 
system 


See  reference  on  page  53  . 


t>v  - 
- 
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THE  ELEMENTS  OF  ORIENTATION 


In  the  Right  Ascension- 
Declination  System  (denoted 
by  *) 


d#  ■  d 
A  X* 

A  y*  sAy 

'  V*  =  90  -  6 

A*  «  360  -  RA 

K* = K+A< 


In  the  Azimuth-Elevation 
System,  j i  latitude  of  station, 
<k  =■  sidereal  time  of  exposure 


d  =  d* 

2  Ax* 

Ay  g  Ay* 

I 

I 

v  »  90  -  £ 

RA  =  0  -  t  (t  =  hour  angle) 

cos  fi  sin  A 

sin Af(  *  - 

cos  6 


sin  6  *  gin  fi  cos  v  -  cos  fi  sin  v  cos  A 

sin  v  sin  A 

sin  t  =  — - 

cos  6 


V.  THE  ACCURACY  OF  THE  METHOD* 

The  accuracy  of  a  measuring  method  is  determined  by  the  propagation 
of  the  systematic  instrumental  errors  and  the  random  errors  of  the  obser¬ 
vations.  The  basic  requirement  of  any  measuring  method  is  that  the  sys¬ 
tematic  errors  *ould  be  sufficiently  small  so  that  their  influence  on  the 
result  can  be  neglected  or  the  systematic  errors  must  be  known  accurately 
enough  in  order  to  apply  corresponding  corrections.  The  systematic  errors 
of  an  instrument  or  measuring  method  may  be  analyzed  by  discussing  critically 
the  assumptions  made  in  designing  the  instrument  and  in  deriving  the  mathe¬ 
matical  analysis  of  the  measuring  method. 


A  detailed  study  of  the  accuracy  of  intersection  photogrammetry  will 
be  published  in  a  separate  report.  See  reference  note  on  page  3U» 
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The  following  criteria  characterize  the  systematic  errors  which  we 
must  deal  with  in  the  solution  of  our  problems 

1)  The  photograph  produced  by  the  camera  should  be  an  exact 
central  projection.  To  achieve  this,  the  measured  plate  coordinates  must 
be  corrected  for  distortion.  The  determination  of  this  distortion  is  a 
part  of  the  camera  calibration.  It  is  not  discussed  in  this  report.  Modern 
measuring  lenses  are  practically  distortion  free,  e.  g.,  the  "Wild  "Aviotar 
lens'1,  with  a  field  angle  of  60°,  has  for  all  apertures  from  lsU.2-to  I2I6 

a  maximum  distortion  of  less'  than  h  microns.  In  general,  the  distortion 
curve  for  a  certain  lens  may  be  expressed  accuract^ly  enough  as  a  function 
of  the  radial  distance  from  the  principal  point  of-  the  plate  and  corres¬ 
ponding  corrections  may  be  applied  directly  to  the  comparator  measurements 
of  the  plate  coordinates. 

2)  The  comparator  measurements  of  the  plate  coordinates  should 
be  affected  only  by  random  errors  in  setting  and  reading.  Hence,  it  is 
assumed  that  the  comparator  is  adjusted  and  calibrated.  The  procedure  will 
not  be  discussed  in  this  report.  Especially  careful  attention  must  be  paid 
to  the  perpendicularity  of  the  mechanical  axes  of  the  comparator  and  to 
the  consistency  of  the  scales. 

3)  The  formulas  for  astronomical  and  terrestrial  refraction 
must  be  adequate.  The  directions  of  the  control  points  are  affected  by 
lateral  and  vertical  refraction.  The  lateral  refraction  is  neglected. 

The  astronomical  refraction  is,  hence,  determined  as  a  function  of  the 
elevation  angle  and  the  temperature  and  pressure  at  the  time  and  location 
of  exposure.  Local  variations  of  the  refraction  coefficient,  if  known, 
may  be  taken  into  account.  Results  obtained  recently  from  a  wide  variety 
of  sources  in  extensive  research  on  precision  trigonometric  leveling  nets 
indicate  that  the  refraction  anomolies  are  less  than  assumed  heretofore, 

U)  The  photographic  emulsion  must  represent  a  plane  with  suf= 
ficient  accuracy  and  no  irregular  shifts  of  emulsion  must  take  place. 

?)  It  Is  assumed  that  the  control  points  are  essentially  free 
of  errors.  The  reduction  coefficients  of  the  RA  and  £  values  for  date 
and  time  of  exposure  are  taken  from  the  American  Ephemeris.* 

6)  The  time  interval  between  the  orientation  exposures  of  the 
different  stations  must  be  known  with  sufficient  accuracy.  (+_  0.01  sec  = 

+_  0.1?  sec  of  arc.) 

7)  The  target  points  must  be  recorded  simultaneously. 


* 


Compare  remark  on  page  ?0  . 
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8)  The  relative  spatial  situation  of  the  camera  stations  must 
be  known  with  sufficient  accuracy  from  a  local  geodetic  net. 

9)  The  cameras  must  not  change  their  interior  or  exterior 
orientation  during  the  time  interval  between  the. individual  orientation 
exposure  and  the  final  target  point  registration. 


The.  triangulation  results  of  the  outlined  photogrammetric  method 
will  be  the  less  affected  by  systematic  errors  the  better  the  above- 
mentioned  conditions  are  satisfied.  If  all  steps  in  the  arrangement  of 
a  photogrammetric  measurement  are  done  with  the  necessary  care  the  results 
will  be  affected  only  by  the  random  errors  of  the  plate  measurements . 

The  accuracy  of  the  plate  coordinate  measurements  depends  on  the  image 
quality.  Star  images  and  additional  recorded  target  points  are  in  general 
measurable  within  a  few  microns.  The  influence  of  such  a  random  error 
(ra)  on  the  corresponding  ray  may  be  expressed  with  sufficient  accuracy  by 
the  relation  of  m  to  the  principal  distance  (d)„  Therefore,  assuming  a 
mean  error  of  a  reading  of  +  3  p,  and  d*-"  300  mm.  we  have  to  expect  a  rela¬ 
tive  angular  deviation 


300,000 


2» 


The  orientation  of  the  plate  will  be  obtained  more  accurately,  when' 
more  stars  are  carried  in  the  least  squares  adjustment.  However,  it  should 
be  realized  that  even  with  an  accurately  determined  plate  orientation, 
the  single  spatial  ray  to  any  target  point  will  be  affected  by  the  entire 
amount  of  the  error  of  the  target  image.  Hence,  it  is,  sufficient  to  carry 
in  the  least  squares  adjustment  of  the  orientation  so  many  stars  as  are 
necessary  to  reduce  the  mean  error  of  the  plate  orientation  to  an  insig¬ 
nificant  fraction  of  the  angular  error  of  the  target  images. 


In  general,  it  may  be  assumed  that  with  a  mean  error  of  +_  2-3  p. 
and  a  focal  length  of  300  mm,  the  individual  direction  in  space  will  be 
obtained  from  a  10-star  orientation  adjustment  with  an  angular  accuracy 
of  1  to  2  seconds  of  arc*  Hence,  the  accuracy  of  the  individual  photo¬ 
grammetric  camera  is  1; 100 '000  to  Is 200  000.  The  propagation  of  this 
error  during  the  triangulation  procedure  depends  on  the  geometry  of  the 
particular  configuration. 


VI.  NUMERICAL  EXAMPLES 


The  first  three  examples  demonstrate  the  validity  of  the  geometry 
of  the  solution  expressed  by  the  formulas  10a,  10*5,  and  14  in  connection 
with  19.  For  each  of  the  three  solutions,  a  rigorous  three-star  computa¬ 
tion  was  computed.  The  least  squares  adjustment  of  an  over-determined 
solution  is  shown  by  three  further  examples.  For  all  computations,  the 
coordinate  system  was  oriented  with  respect  to  the  horizon  and  Zenith 
of  the  station. 
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1.  Example  based  on  the  formula  10a  of  the  projection  method. 

We  use  the  results  of  a  former . calibration  as  the  following  approximations  s 


(length  unit  is  1  decimeter  =0.1  meter) 


=  +.00192 

Sin  a*  ■  +  .00156110 

0 

A  -  -.72967230 

=  -.00188 

cos  tc ^  =  +.99999878 

B  -  -.  6303*014; 

=  +3.01113 

sin  vQ  »  +.3U0999l;5 

C  -  -.265038# 

cos  v  ■  + . 940063 1;9 

0 

A*  =  -.59270391; 

=  +0°  51  22" 

B*  -  +.77631620 

-  19°  56'  16" 

sih  =  -.62920358 

C«  «  -.211;55807 

=  2180  59’  29" 

cos  Aq  ■  -.7772l;Q53 

D  »  +  . 3*4)99903 

E  -  +.00053233 

F  =  -,9i;00  63 1;9 

The  standard  coordinates  of  the  three  stars  computed  from  astronomical 
data  with  formula  (2)  ares 

^  c 

3  +.16900891  +.OU650153 

10  +.15713779  +.38332881 

18  +.l;8l27li91  +.396l327i<. 


I  £'  JL-bx.  j£'-Ay  s  t 

,  o  0 

3  +0.21350  -0.^7731  +0.21158  -0.575U3  -0.5115583 9  -0.1U07U96U 

10  -o.56ll;5  +0.00056  -0.56337  +0.00214;  -0.IM729U8  -1.08U86171 

18  -0.01032  +0.63807  -0.01221;  +0.63995  -l.25730l;Olt  -1.03I486525 


£  -  SL 

o  u 


T7  =  -  <£  _<£  A  77=77  -n 

'o  u  c  o  '  'c  '0\ 


3  -?3 .02676182  +.16901178  +  »Oi|650172  -.00000287  -.00000019 

10  -2.83012121;  +  .1571201*9  +.38332691;  -.00000370  +.00000187 

18  -2.612U3757  +.U8127621  +.3961301U  -.00000130  +.00000260 


a 


a1- 


3  +.00006172.-  -.00000112  -17.29801637  +3.7670319*; 

10  +.00000965  +.00001190  -  2.17656068  +2; . 8903U92I; 

18  +.00000328  +.ooooo51;0  -  2.71937570  +1.820U0370 
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3 

10 

18 


-3.76092138 

+0.1*08161*1*3 

+0.1*0052769 


-0.62151661* 

-0.91211*623 

+0.38791*159 


+1*.  1*7911206 
+0.58111286 
+0.609351*35 


+1.517503H 

+1.71*513312 

+0.61727737 


+5.59365660 

+0.72198821* 

+0.863671*78 


b} 


c 


c1 


3 

10 

18 


-1.189,62171* 

-1.62661*71*0 

-0.57881*1*82 


-04751*22l*2l* 
-0.10813, 877 
+0.18107826 


+0.33396677 

-0.32782533 

-0.12553003 


[a  +  a«]  /2n  =  -1.952691*61* 
[P  +  P'3  /2n  =  -0.68299176 
[a  +  a'J  /2n  -  +1.59158211* 
Lb.+  b']  /2n  =  +0.63070091* 
[c  +  c']  /2n  =  -0.1331*1*556 
|L  +  L'3  /2n  =  +0.000015H* 


A 

Uv) 


B 


C  D  E 

(Mx)  (A  Ay)  (Ad) 


L 


3-  -15.31*532173 
10  -  0.22386601* 
18  -  0.76668106 
3  +  5.71972658 
10  +  6.81*301*388 
18  +  3.77309831* 


-3.07792962 
+1.09115619 
+1.08351 91*5 
+0.0611*7512 
-0.229151*1*7 
+1.07093335 


+2 .88752992 
-1. 0101*6988 
-0.98222779 
-0.071*07903 
+0.15355098 
-0.971*301*77 


+1*. 96295566 
+0. O9128730 
+0.23297381* 
-1.82032268 
-2.25731*831* 
-1.20951*576 


-0.62077868 

+0.02530679 

+0.311*52382 

+0.1*671*12335 

-0.191*37977 

+0.00791553 


-0.00001*658 

+0.0000051*9 

+0.00001186 

+0.00001626 

+0.00000321* 

+0.00000971* 


The  normal  equations  are:* 


CAv) 


(A*)  (AAx) 


(AAy) 


+329.8956076  +1*8.9810773  -1*6.3799138  -106.7798210 

...+13.01*11*762  -12.1376130  -  15.8135538 

...+11.301981*0  +  11*. 9763025 

...♦  31*.  5657362 

(Ad)  (L) 


+10.6521*313 
+  2.3608737 
-  2.1992076 
.  -  3.1*2691*37 
...♦  0 .71*12523 


+  .0008563871*18 
+.000172898966 
-.000161891*590 
-.000276603 1*69 
+  .0000398321*89 


0 

0 

0 

0 

0 


« 


Weighting  factors  are  omitted  because  this  is  a  unique  solution. 


Reduction  of  the  normal 
equations  gives? 

AA  =  +.00000793  -  +1.6" 
Av  -  +.00000584  =  +1.2" 
-.00000646  -  -1.3" 
Ad  =  -.00002166 
AAx  «  -.00000144 
A  Ay  B  +.00002157 


The  final  orientation 


elements  are: 

A  =  218°  59*  30.6" 
v  =  19°  56'  17.2" 
X=  +0°  5'  20.7" 
d  =  .30111083 
Ax  =  +.00019186 


(Length  unit  1  meter) 
sin*®  +  .00155480 
cosx"=  +  .99999879 

sin  v=  +.341001*92 
cos  v=  +.94oo6l5l 


Ay  =  -.00018584  sin  A=  -.62920961 

cos  A-  -.77723565 


c 


The  check  is 

obtained  with  formula  10. 

X-Ax 

£  -£ 

^c  ^  7} 

V’ 

in  8th 

in  8th 

decimal  place 

decimal  place 

3 

+0.02115814 

-0.05754516  +.16900894 

.  -3  +.04650148 

+5 

10 

-0.05633686 

+0.00024184  +.15713778 

+1  +.38332876 

+5 

18 

-0.00122 386 

+0.06399284  +.48127491 

0  +.39613266 

+8 

2,  Example  based  on  the  formula  10°  of  the  projection  method. 

We  use  the  results  of  a  former  calibration  as  the  following  approximations* 
(length  unit  is  1  decimeter  =0.1  meter) 


Ax  ■  +.00192 
o 

AyQ  =  -.00188 

d  =  +3.01113 
•o 

Ko  =  52°  38’  20.02" 
coo  =  12°  23'  22.57" 
aQ  =  1*°  44'  42.37" 


sin  XQ +.75482675 

cos  k  =  + .60683642 
0 

sin-w  =  +.21455807 
COS  05q  ■  +.97671123 

sin  aQ  =  +.27135819 
cos  aQ  =  +,96247843 


A  =  +.72967232 
B  =  +.63034344 
C  =  +.26503859 

A'=  +.59270395 
B’=  -.77631621 
C»-  +.21455807 

D  =  -.34099903 
E  =  -.00053234 
F  =  +.94006349 


The  standard  coordinates  of  the  three  stars  computed  from  astronomical 
data  with  formula  (2)- ares 

c  ^  c. 

3  +.16900891  +.04650153 

10  +.15713779  + .38332881 

18  +.48127491  +.39613274 
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Measured  plate 

coordinates 

l 

l'  i-dx 

0 

0 

* 

I 

s  t 

3 

+0.21350  -0.57731  +0.21158 

-0.5751*3 

+0.51155839  +0.11*071*961* 

10 

-0.5611*5  +0.00056  =0.56337 

+0.0021*1* 

+0  441*7291*8  +1 .081*86171 

18 

-0.01032  +0.63807  =0.01221* 

+0.63995 

+1.257301*01*  +1.031*86525 

u 

O  u 

„  t 

V  *  — 

'0  u 

*v-\-va 

3 

+3*02676182 

+.16901178  + 

.01*650172 

=.00000287  -,00000019 

10 

+2.83012121* 

+.15711*11*9  + 

.38332691* 

-.00000370  +.00000187 

18 

+2.6121*3757 

+  .1*8127621  + 

,39613011* 

-.00000130  +.OOOOQ260 

+L«- 

0 

0 

3 

-.00000279 

=  .000021*18 

10 

-.00000361 

+  .00003101* 

18 

-.00000106 

+  .00001361* 

a 

a* 

0 

ps  a 

3 

-.001*91382 

+128.5632671*  ■ 

=.17003235 

+19 .3651*5308  -.20250197 

10 

-.01*1*93296 

+  19.031*0551*  J 

►. 15268999 

+  2.3795511*9  -.21738815 

18 

+.06170817 

+  6.5256801  H 

►  .12882238 

=  0.97932866  =,19598763 

a" 

b 

b' 

c  c' 

3 

+3^63326183 

=.25289086  =25.58231098 

+  .031*09870  +7 .181811*25 

10 

+  1^5259701* 

=  .27008811*  - 

1*.  21*31*9876 

+  .0l*04&$£  -0 .85521071 

18 

+1.55826911* 

=  .27778511*  - 

1.1*6121*911* 

-.05821*050  -0.31689090 

+  a'J  /2ri  =  +25.68911*1*05 
[P  +  0*3  /2n  =  +2,1*7952599 
[a  +  a»l  /2n  *  +6.351*70838 
[b  +  b»3  /2n  -  =5.31*797050 
Jc  ♦  cQ  /2n  -  +1.001*3371*0 
[L  +  L»]  /2n  =  +.00000217 
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L 


A  co 5  •  10  A  co 

ABODE 
(A<b)  (Art)  (lAx)  (AAy)  (Ad) 

3  -2.5691*0578  -2.61*955831*  -6.55721035  +5.09507961*  -0.97023870  +  .000001*96 

10  -2.573U0770  -2.32683600  -6.57209653  +5.07788236  -0.96388381*  +.00000578 

18  -2.56271*358  -2.35070361  -6.55069601  +5.07018536  -1.06257790  +.00000323 

3  +10.2871*1231;  +10.8859270 9  +26.2785531*5  -20.231*31*01*8  +  6.1771*7685  +.00002635 
10  -  .66550886  -0.099971*50  -1.80211131*  +1.101*1*7171*  -1.85951*811  -.00002887 

18  -1.91631*639  -3.1*58851*65  -1*. 7961*3921*  +3.88672136  -1.32122830  -.0000111*7 

The  normal  equations  are:* 

(Ago)  (A^C)  (A  Ax) 

+129.7380661  +137.5028570  +331.27781*73 

...+11*8.1*372116  +350.9011*152 

...+81*5.916881*7 


(Ad) 

+71.01621*61  +.000276370838  =  0 

+79.311*7160  +.000295219811  =  0 

+191.6808813  +.000707812766  =  0 

-11*7.1*116336  -.00053861*3211  =  0 

...+  1*6.361*2906  +.000217800398  =  0 


Weighting  factors  are  omitted  because  this  is  a  unique  solution. , 


(A  Ay) 

+255.1*91*7177 

-271.0571917 

-652.3570867 

,..+503.2065008 


Reduction  of  the  normal 
equations  gives: 

Aa  =  +.00000295  =■  +0.61" 
Aco  =  +.00000577  =  +1.19" 
tJC  =  -.00000038  -  -0.08" 
Ad  =  -.00002163 
A  Ax  =  -.00000187 
AAy  =  +.0000211*0 


The  final  orientation 


elements  are: 

a  =  15°  1*1*'  1*2.98" 
co  =  12°  23'  23.76" 
K  =  52°  38*  19.91*" 
d  -  .30111081* 

Ax  «  +.00019181 
Ay  =  -.00018586 


(Length  unit  1  meter) 

or  transformed  with 
formula  (ll) 

A  =  218°  59*  30.6" 

v  *  19°  56'  17.2" 

(*)  =  +  0°  5e  22.7" 

sin  co  =  +.211*56372 
cos  co  =  +.97670999 


sin  a  *  +.27136103 
cos  a  =  +.9621*7763 


The  check  is  obtained  with  formula  lO*3. 


sin  tC  -  +  .791*82651 
cos  K  -  +.60683673 


61; 


V-Ay 


&€  * 

-din 

c 


A  7  “ 

?c~7in 


3 

+0.02115819 

“0 .05751+511+ 

8th  dec.  place 

+.16900893  -2 

+  .ol+65oi55 

8th  dec.  place 

-2 

10 

-0.05633681 

+0.00021+186 

+.15713778  +1 

+ .38332879 

+2 

18 

-0.00122381 

+0.06399286 

+  .1+81271+89  +2 

+  .  39613276 

-2 

3.  Example  based  on  the  formulas  (lit)  and  (19)  of  the  analytical  method. 
We  use  the  results  of  a  former  calibration  as  the  following  approximations: 
(length  unit  is  1  meter) 

a1°  =  +2. 2281+2997 
b1°  -  +2,57822963 


Cl°  =  +0,28198337 
a2°  =  -2.71+202238 

b2°  =  +2.091+37269 
c2°  =,'+0.229131+71+ 


And  the  auxiliaries  are: 

a’  =  +.06119181+ 
b'  =  -1.0106391+8 
c'  =  +.26700761+ 
d'  =  +.231611+68 


a  0  =  +0.00211+801 
0 


bQ°  =  -1.201+51+1+09 


These  approximations  satisfy  the  condition  equations  (19) . 


l)  a  °b  0  +  a,  Ub1 
00  1  1 


°^  °  +  a2°b2°  =  +.O0t)OOOOl6 


2  2  2  2  2  2 
00  0  ,  o  ,0  ,  o 

2)  ao  +al  *  a2  -  bo  "bl  -b2 


+  .000000023 


The  measured  plate  coordinates  1  and  1?  are  corrected  for  distortion  and 
comparator  constants . 

(1)  (2)  (3) 

a.^%  t-b^/’+c^0  a2CX+b2°£>+c2  aQ°£+boCjZ'+l 


l  V 

(meters) 

3  +.021350  -.057731 
10  -.05611+5  +.000056 
18  00103 2  +.063807 


+  .18071658  +.01+968233  +1.06958539 

+.15701255  +.38320287  +0.99981195 

+  .1+1+1+19273  +.36560015  +0.9231391+1+ 
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Computed  from  formula  (2) 


(b)  (5) 

“•c  “9 c  X  (3)  X  (3) 


(meters) 


3 

.  + .16900891 

+. 0b 6501 53  + .18076946  . 

+.0b973736 

10 

+.15713779 

+.38332881  +. 1571082b 

+.38325673 

18 

+.b8l27b91 

+. 3961327b  +.bbb28385 

+  .36568576 

-L=(l) -(b) 

-L»=(2)-(5)  1 v 

I*? 

(meters) 

3 

-.00005288 

-.00005503  +.00099281 

-.00268458 

10 

-.00009569 

-.00005386  -.02192200 

+ ,000021b7 

18 

-.00009112 

-.00008561  -,000b0881 

+.02 52760b 

a 

P  r 

5 

3 

+.05965167 

+.0l80bbb9  .-.01127907 

-,0106025b 

10 

-„003b9221 

-.05673893  -.00615567 

+ .01301016 

18 

-.06b5b902 

+.00286lb9  +.01b99976 

-.013b3273 

C<0/3 

=  -.00279652 

Cpl/3 

=  -,0ll9bb32 

Zrl  /3 

=  -.00081166 

r  ©3  /3 

-  -.0036750b 

-L 


+.00002702 
-.00001?79 
-.00001122 
+.00000980 
+.00001097 
-.00002078 

*  The  decimal  point  in  the  coefficient  A...D  was  moved  for  the  convenience 
of  the  numerical  computations. 


[-1-  1  /3 

[11/3 

r^'D/3 

Uvi  n 

[i'  /3 

C-L']/3 


=  -.00007990 
«  -.0119U233 
=  +.0020bb00 
=  -.00697933 
=  +. 0075376b 
=  -.0000 6b 83 


The  observation  equations  ares* 


100  ^a2^b2 


(^a2) 

3  +.06.2bb8l9 
10  -.00.069569 
18  -.06.175250 
3  +  .03 .32^233* 
10  -,0b.b20267 
18  +.01.091033 


E100Ab2)Aao 

B 


Too  Aao)*bo  =100^bo 


(Ab2) 

+.02.998881 

-.0b.b79U6l 

+.0l.b8058l 

-.05.977500 

-.00.193800 

+.06.176300 


C 

(Mq) 

-.01.0b67bl 

-.00.53bb0l 

+.Ol.58llb2 

-.00.79721b 

+.01.b5b267 

-.00.657052 


(Abo) 

-.00.692750 

+.01.668520 

-.00.975769 

+.01.022222 

+.00.751617 

-.01.7738b0 
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The  normal  equations  ares  * 


da*2 

Ab«2 

Aa» 

0 

Ab* 

0 

+108.91*92225 

-2.3870512 

-26.0627025 

-0.2710087 

+  .000200581*298  =  0 

• 

..+105.1676876 

+2.0138737 

-28.2117716 

-;oooo5395£?i£  =  0 

...+7,0631*212 

-0.  2657310 

-.000015791011  =  0 

... +8. 9723610 

+.000021007502  =■  0 

*  Weighting  factors  can  be  neglected  because  this  is  a  unique  solution. 
The  solution  of  the  normal  equations  is; 


Aa'  = 
0 

Ab'  = 
0 

-.00001*031*01* 

-.000001*2328 

or 

A  a  =  -  .00l*03l*0l* 
0 

Ab  =  -.0001*2328 
0 

Aa*2  - 

Ab'2  = 

-.000011501*2 

-.0000001111 

• 

Aa2  =  -.0011501*2 
Ab2  =  -i 00001111 

And  with  the  formulas  (1*1*)  and  (50)  s 


=  -.00106071 

Abx  ~  +,00011462 

Acx  =  +.00007172 
Ac2  =>  +.00007608 


The  plate  constants  are  now; 

a1  =  +2.22736926 

\  =  +2.57831425 
^  =  +0.2820550 9 
a2  =  -2.71*317280 
b2  =  +2.091*36158 
c2  =  +0.22921082 
aQ  -  -0.00188603 
bQ  -  -1.201*96737 


The  condition  equations  reduce  to 

1)  + .00000161 
2)  +.00001853 
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If  the  solution  is  repeated  to  eliminate  the  influence  of  second  order 
terms,,  we  obtain  the  following  new  L  and  L’  terms t 

-L  -L» 

3  +.00000075  +  .00000000 

10  +.00000089  +  .0000001*6 

18  +.00000132  +  .00000063 

and  the  new  absolute  terms  in  the  normal  equations  ares 

TA(L)]  -  ~.OOOOOli87717U 

•  [b(L)3  =  +.00000U01633U 

[C(L)1  -  +.000001078512 

[D(L)]  =  -.000001090959 

The  corresponding  reduction  gives* 

Aa1  =  +.00000127 
i^b1  =*  -.00000399 
Ac1  =  -.00000098 
-  +.00000856 

Ab2  -  -.00000312 

Ac2  -  -.00000035 

Aa  «  +.00001732 
0 

Ab  =  +.00000313 
0 

And  the  final  plate  constants  ares 

a-L  =  +2.22737161  a2  =  -2 .714316815  aQ  =  -0.00187329 

\  -  +2.57831^83  b2  -  +2. 09^6279  bQ  -  -1.201*9651*8 

c1  -  +0.28205510  c2  =  +0.22921080 

The  condition  equations  reduce  now  tos 

1)  +.000000025 

2)  -.000000069 
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The  final  check  with  equations  (lU)  gives; 


l 

V 

< 

A<t  * 

A  7}  * 

3 

+.021350 

-.057731 

+.16900890 

+1 

+.OU650153 

0 

10 

-.0*6lU5 

+ .000056 

+.15713779 

0 

+.38332882 

-1 

18 

-.001032 

+.063807 

+  .U8127 h92 

-1 

+  .39613273 

+1 

The  differences  A£  and  A  7)  are 

in  units 

of  the  8th 

decimal 

and  show  the  difference  against  the  originally  computed  standard 
coordinates » 

The  computations  of  the  orientation  elements  with  formula  (16)  gives; 

tan  PC  =  +  ,001**l;6l; 

|<  =  +0°  *>  20^ 

sin  PC  -  +  .00155U61; 
cos  /c  *>  +.99999079 

7f  .  =  -.13661U81 

-  -.110*9631 

u 

cot  (A+<)  =  +1.2313373* 

(A+x)  -  219°  U'  51.31" 

sin  (A+X)  -  -.63010.732 
cos  (A+tf)  =  -.776256^0 

A  =  2l9u  *9?  30. 6U" 

sin  A  =  -.62920976 
cos  A  »  -.77723553) 

h‘  =  +2.93250068 

xc'  -  -.00019211; 
y  »  -  +.05311195 

d'  =  +.88301011; 
sin  v  =  +.3U100U96 
lv  °  19°  17.21" 

v/2  »  9°  58'  8.60" 

tan  v/2  -  +.17577016 
1  d  =  .30111081; 
y  '  »  +.05292630 
y0*  -  yc.. »  -  -.00018555 

Ax  =  +.0001918* 

Ay  -  -.00018585 
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1*.  Example.  Least  squares  solution  based  on  formula  10a  of  the 
pro jsctii.on  system. 


We  use  the  results  of  . 

a  former  calibration  as 

the  following  approxi 

mations .  # 

(length  unit  is  1 

decimeter  =  0.1  meter) 

Axq  =  -.00050 

sin  (CQ  =  -.00058178 

A  =  -.72953695 

AyQ  =  -.00150 

cos -  +.99999983 

B  «=  -.63070287 

d  -  3.01110 

O 

sin  Aq  =  -.63112719 

C  »  -.261*55573 

cos  Ao  =  -.77567936 

A*  =  -.59283365 

OvJ 

0 

0 

1 

D 

*C 

B«  =  +.776021*39 

A  =  219°  8' 

0 

sin  vQ  *  +.3U106326 

C'  =  -.215251*30 

V  =  19°  56*  30" 

0 

cos  v  =  +.91*001*035 

0 

D  *  +.31*106320 

E  =  -.0001981*2 

F  =  -.91*001*035 

The  standard  coordinates  of  the  four  stars  computed  from  astronomical 
data  with  formula  (2)  are; 

£  V 

3  + . 16900891  + . 0^650153 

10  +.15713779  +.38332881 

17  +  .51*637688  +.15537271 

18  +  .1*81271*91  +.3961327^ 

Measured  plate 
coordinates 


8, 

V 

s 

t 

3 

+0.21350 

-0.57731  +0.211*00 

-0.57581 

-0.5111*995 

-0.11*07231* 

10 

-o.56ll*5 

+0.00056  -0.56095 

+0.00206 

-0.1*1*1*3139 

-1.081*681*3 

17 

+0.60320 

+0.1*0158  +0.60370 

+0.1*0308 

-1.1*711*209 

-0.1*186257 

18 

-0.01032 

+0.63807  -0-,  00982 

♦0.63957 

-1.2570002 

-1.031*9311* 

A  7) 

u 

^0 

Vo 

c  0 

"7  -7 

• c  '0 

3 

-3.0269856 

+.16897983 

+  .01*61*8962 

+.00002908 

+  .00001191 

10 

-2.82971*16 

+.15701571 

+.38331567 

+ .00012208 

+  .00001311* 

17 

-2.6931995 

+  .51*631*679 

+  .1551*3806 

+  .00003009 

-.00006535 

18 

-2.6121*198 

+  .1*8116317 

+.39615815 

+.0001117U 

-.000025U1 
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+l=t„4.£  +Lia,lA$ 

% 

3  -.00062552  +  .0000701*8 

10  -.0003181*8  +.00008368 

17  -.00019358  -.00011961 

18  -.00028206  -.00005281 

a  a*  p  a 

3  -17.26807537  +3.7790591*3  -3,77762362  -0.61891*165  +1*. 1*8162381 

10  -2.17259130  +1*.  9071*3739  +0.1*0617633  -0.90991300  +0.5811*331*3 

17  -6.821*66631  +1.30365721  -0.71360681*  +0.1*7757286  +1.50631*1*36 

18  -2.711*99269  +1.821*81662  +0.398181*03  +0.3891991*1  +0.60932290 

a*  b  b'  c  c* 

3  +1.51717371*  +5.5937371*3  -1.19001006  -0.751171*79  +0.33539071 

10  +1.71*673907  +0.72195129  -1.62850750  -0.10782366  -0.32652118 

17  +0.5271*1892  +2.18781531*  -0.1*3892801  +0.591*88058  +0.01*698605 

18  +0.6171*21*72  +0.8631*8139  -0.57  9U.576  +0.l8ll*19Ql  -0.12502015 


P  P' 


3 

0.0195 

0.2690 

10 

1.1638 

0.1753 

17 

0.11*17 

2.1239 

i8 

0.8951* 

1.3955 

La+a*]  /2n  =  -2.11*566938 
[p+p»j  /2n  =  -0.51*361906 
[a+a  /2n  =  +1.1*1*81*3512 
(b+b']  / 2n  =  +0.69130302 
[c+cO  /2n  =  -0.01898283 
[L+L»l  /2n  -0.00017971* 


ABC 
(Av)  (AAx) 

3  -15.1221*0599  -3.231*001*56  +3.03318869 

10  -0.02692192  +0.91*979539  -0.86700169 

17  -1*.67899693  -0.16998778  +0.05790921* 

18  -0.56932331  +0.91*180309  -0.83911222: 

3  +5.921*72881  -0.07532259  +0.06873862 

10  +7.05310677  -0.36629391*  +0.29830395 

17  +3.1*1*932659  +1.02119192  -0.92101620 

18  +3.9701*8600  +0.9328181*7  -0.8310101*0 


D  E 

U Ay)  (A  a) 

+1*. 9021*31*1*1  -0.73219196  +.0001*1*578 
+0.03061*827  -0.08881*083  +.00013871* 
+1.1*9651232  +0.6138631*1  +.00001381* 
+0.17217837  +0.2001*0261*  +.00010232 
-1.88131308  +0.351*37351*  -.00025022 
-2.31981052  -0.30753835  -.0002631*2 
-1.13023103  +0.06596888  -.00006013 
-1.2701*1878  -0.10603732  -.00012693 


The  normal  equations  ares 

(&v)  (4*)  (AAx)  (AAy)  (Ad) 

+73.2853591  +12.6331*215  -11.3511*01*3  -23 .7121*1*01*  +0.2101651*  -.0020651*3059?  «  0 

+5.50631*56  rl*. 9585661  -1*.  081*1927  +0.1197973  -.0000621*63830  =  0 

+1*. 1*673729  +3.6700301  -0.111*801*1  +.000056030968  »  0 

+7.671*5081*  +0.0632622  +  .0006693981*1*8  -  0 

+0.181*2911*  -.6000001*38911  »  0 

+  .000000091*81*7  -  0 
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Reduction  of  the  normal  equations 
gives? 

A  A  -  -.00007911  -  -16. 3 “ 

A.v  =  +.00000627  -  +1.3“ 

A>c=  -.00011*390  «  -29.7" 

Ad  =  +.00010869 
A  Ax  =  -.0000^625 
A  Ay  =  -.0001181*1* 

£  +  .0000000081*10 


The  final  orientation  elements  are: 
(Length  unit  1  meter) 

A  =  219°  7'  1*3.7" 

v  =  19°  56'  31.3“ 

*=  - 0°  2'  29.7" 

d  =  .30112087 

AX  =  -.OOOC5562 

Ay  =  -.00016181* 


The  computation  of  the  residuals  with  formulas 


(2U)  and  (26): 


P 

P' 

V 

v' 

3 

-all.  M3 

+55.218 

+21.11 

-19.1*9 

10 

+37.377 

+58.01*1 

+1*6.55 

+11*.  29 

17 

-10U.80I* 

+7.375 

-16.1*6 

-36.57 

18 

+11.811* 

-50.596 

-1*1.59 

+1*3.03 

Ivppl  - 

81*10. 

I Vv}  = 

8386. 

The  final  check  is  obtained  with  formula  (20) 


x-  y=  V  +  v  1 

<£ 

V 

8th  dec 

3  +.02135211  -.05773295 

+ .16900889 

+2 

10  -.05611*031*  +.0000571*3 

+.15713778 

+1 

17  +.06031835  +.01*0151*31* 

+  .51*637681* 

+1* 

18  -.00103616  +.06381130 

+  .1*81271*90 

+1 

The  mean  error  of  an  observed 
1  or  l1  value  is: 


of ter  the  residuals  are  applied. 
4»  » 

’  %-* 

8th  dec. 

+  .01*650156  -3 

+.38332882  -1 

+.15537275  -U 

+.39613276  -2 


A  certain  systematic  deviation  in  the  individual  parameter  of  the  solution 
must  be  expected  due  to  the  very  unfavorable  weight  p  of  point  3.  This 
small  number  is  caused  by  the  procedure  which  leads  to  the  elimination  of 
A A,  In  the  next  example  this  reduction  step  was  not  carried  out  arid 
consequently  the  normal  equation  system  has  six  unknowns,  but  more  nearly 
even  weights.  A  comparison  of  the  results  of  these  least  squares  adjust¬ 
ments  at  the  end  of  this  paragraph  shows  the  influence  on  the  numerical 
values . 
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5.  Example,  Least  squares  adjustment  based  on  the  formula  10  of 
the  projection  method. 

We  use  the  results  of  a  former  calibration  as  the  following  approxi¬ 
mations  s 

(length  unit  is  1  decimeter  *?  0,1  meter) 

Axq  -  -.00050  sin  -  +.79465232  A  =  +.72953695 

Ayo  =  -.00150  cos  k  =  +.60706445  B  *  +.63070288 

c  -  +.26455573 

d  =3,01110  sin  co  =+.21525430 

o  o 


cos  co  =  +.97655803  A!=  +.59283366 

o 


*o  -  52° 

375 

20.84" 

B»-  -.77602439 

ti)  -  12° 

o 

2$* 

49.61" 

sin  a  =  +.27090631 
o 

C*=  +.21$2$430 

«o  * 

43' 

$.54" 

cos  a  o  ♦ ,96260572 

0 

D  =  -.34106320 
E  =  +.00019842 
F  =  +.9400403$ 


The  standard  coordinates  of  the  four  stars  computed  from  astronomical 
data  with  formula  (2)  are5 


*c 

*c 

3  +..169O0891 

+ .04650153 

10  +.15713779 

+.38332881 

17  +.54637688 

+ .15037271 

18  +  .48127491 

Measured  plate 
coordinates 

+.39613274 

JL  V 

h-bx 

0 

3 

+0.21350  -0.57731 

+0.21400 

10 

-0.56145  +0.00056 

-0.56095 

17 

+0.60320  +o.4oi58 

+0.60370 

18 

-0.01032  +0.63807 

-0.00982 

t'-4yo  s  t 

-0.57081  +0.511499$  +o. 1407234 
+0.00206  +0.4443139  +1.0846843 
+0.40308  +1.4714209  +0.4l862$7 
+0.639$7  +1.2570002  +1.0349314 


u 

*0 

+L'  =  >?  -h 
'c  '0 

3 

+3.0269856 

+.16897983 

+ .04648962 

+.00002908 

+ .00001191 

10 

+2.8297416 

+.15701571 

+.38331567 

♦ .00012208 

+.00001314 

17 

+2.6931 995 

+.54634679 

+.15543806 

+ .00003009 

-.00006535 

18 

♦2.6124198 

+ .48116317 

+.39615815 

+ .00011174 

-.00002541 

a 

a* 

P 

P‘ 

r 

3 

+1.028551*18 

+.00785581 

-.00503228 

+1.0101*6388 

-.17562028 

10 

+1.021*65393 

+.06018658 

-.01*590669 

+1.11*657879 

+.15569375 

17 

+1.2981*91*81 

+  .081*92309 

+.03963830 

+1.13387202 

-.11092166 

18 

+1.23151800 

+.19061671 

+.07616699 

+1.21*1*0281*3 

+  .15771*386 

r' 

a 

a' 

b 

V 

3 

+  .101*58865 

-.20831*898 

+.25637175 

-.26005071 

-.20108770 

10 

+  .11*287061* 

-.22287255 

+  .271*2651*8 

-.27673521* 

-.25570127 

17 

-.2609201*0 

-.231*11*321* 

+  .28815361* 

-.31*006977 

-.23980691 

18 

“.1872681*3 

-.21*138821* 

+  .29708201* 

-.31*207519 

-.27861*918 

c 

c' 

P 

P' 

3 

+  .031*92183 

+  .056671*26 

9.0061 

9.1*195 

10 

+  .01*133050 

-.05126896 

7.9205 

7.1121 

17 

-.0921*6708 

+.02567068 

5.8661 

7.1151* 

18 

-.07187091* 

-.06015509 

5.7050 

6.0276 

The  normal  equations  ares 

(A  a) 

( A  ca) 

(A*) 

( AAx) 

(a  Ay) 

+36.6831*722 

+3-09791*1*8 

-.1*037976 

-6.5663873 

-10.2379357 

e  ©  c 

.+37.50311*66 

-1.351*6933 

+9 .1603103 

-8.1386151 

•  »  o 

,+1.6279077 

-.31*921*29 

+.2853012 

•  < 

,.+3.7152807 

-.0603955 

• 

..+1*.  2  81*6971 

(A<f) 


-.6217198 

-.002212233932 

*»  0 

-.1919527 

+  .0001*871*21823 

=  0 

+.0152028 

-.0003606231*51 

=  0 

+.0692633 

+ .000590375590 

*  0 

+  .1951*602 

+ .000506051392 

=  0 

...+:i795880 

+  .000011*181*051 

=  0 

. . .+  .00000023901*5 

=  0 

Reduction  of  normal 
equations  gives: 

A  a  =  +  .0000191*5  *  +i*.Ol« 

Ao  =  +.00000132  =  +0.27" 


The  final  orientation 
elements  are: 


a  *  l5° 

co  =  12° 


U3* 

25’ 


(Length  unit  1  meter) 

9.55* 

1*9.88" 


+.0002181*6  =  +1*5. 06« 

K  =  52°  38>  5.90" 

Ad  -  +.00011235 

d  =  .30112121* 

Or  transformed 
with  formula  (ll) 

AA  x=  -.00011082 

Ax  =  -.00006108 

A  =  219°  7'  53.65n 

AAyt  -.00009036 

Ay  =  -.00015901* 

v  -  19°  56’  33.28" 

[p^o)=  +.000000008325 

(*)=  -o°  2’  37.53M 
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The  computation  of  the  residuals  with  formulas  (2l+)  and.  (26)  s 


_/> 

V  V* 

3 

+3.061 

+8.551+ 

+H+.82  -23.61+ 

10 

-13.852 

+7.709 

+1+2.71  +15.66 

17 

+17.271 

+1+.119 

-17.78  -38.51+ 

18 

“3.61+2  ■ 

-21+.652 

-1+3.92  +1+1.61+  The  meajl  error  Df  an  observed 

1  or  1*  value  is: 

£p/¥]‘ 

-  8311. 

tw-:  =  8312. 

m  ■ 

■m 

6.5m- 

The 

final  check 

is  obtained 

with  .formula ( (32)  after 

the  residuals 

are 

applied . 

x-i+v 

y*£.'+vr 

<f  *.-« 

4  lj  » 

V* 

3 

+  .02135H+8 

-.05773336 

+.16900893  8th2dec. 

+  .01+650157 

8th^dec . 

10 

-.05611+073 

+.00005757 

+.15713778  +1 

+.38332882 

-1 

17 

+.06031822 

+  .ol+oi51+i5 

+  .51+637688  0 

+.15537273 

-2 

18 

-.00103639 

+ .06381116 

+  .1+81271+91+  -3 

+  .39613273 

+1 

6.  Examples  Least  squares  adjustment  with. four  stars  based  on 
analytical  method  formulas  (ll+)  and  (19). 


We  use  the  results  of  a  former  calibration  as  the  following  approxi¬ 
mations  s 

(Length  unit  is  one  meter) 


a-j°  =  +2.2261+2997 
b1°  =  +2.57822963 
c1°  -  +0.28198337 
a2°  =  -2.71+202238 
b2°  =  +2.091+37269 
c2°  -  +0.229131+71+ 
aQ°  =  +0.00211+801 
bQ°  =  -1.201+51+1+09 


And  the  auxiliaries  are  : 

a«  -  +0.06119181+ 
b'  -  -1.0106391+8 
c«  =  +0.26700761+ 
d»  =  +0. 231611+68 
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These  satisfy  the  condition  equations  (19) 


1)  +.000000016 
2)  +.000000023 


The  measured  plate  coordinates  1  and  1*  are  corrected  for  distortion  and 
comparator  constants . 

1  (meters)  *  *i°£ 


3  +.021350  -.057731  +.18071658  +  .01*968233  +1.06958539 

10  -.0561U5  +.000056  +.15701255  +.38320287  +0.99981195 

17  +.060320  +  .01*0158  +.51993881  +  .11*781*177  +0*9517571*9 

18  -.001032  +.063807  +  .1*1*1*19273  +.315*0015  +0.9231391*1* 


Computed  from  formula  (2) 

© 

© 

?c 

*cr© 

7c*© 

3 

+.16900891 

+  .01*650153 

+  .1807691*6 

+  .01*973736 

10 

+.15713779 

+.38332881 

+  .15710821* 

+.38325673 

17 

+  .51*637688 

+.15537271 

+.52001829 

♦  .11*787711* 

18 

'  ♦  .  1*81271*91 

+  .39613271* 

♦  .1*1*1*28385 

+ .36568576 

-L=  @  -0 

-L«=0-0 

3 

-.00005288 

-.00005503 

+.00099281 

-.002681*58 

10 

-.00009569 

-.00005386 

-.02152200 

+  .0000211*7 

17 

-.0000791*8 

-.00003537 

+.00937208 

+  .0062391*6 

18 

-.00009112 

-.00008561 

-.0001*0881 

+  .02527601* 

a  p 

r 

6 

p 

3 

+.05965167 

+  .01801*1*1*9 

-.01127907 

-.01060251* 

.0787 

.0823 

10 

-.0031*9221 

-.05673893 

-.00615567 

+.01301016 

.0792 

.0711 

17 

-.036891*17 

+  .0631*1911 

-.00755037 

-.02518991 

.061*8 

.0786 

18 

-.061*51*902 

+  .0028611*9 

+  .011*99976 

-.0131*3273 

.0670 

.0707 

Cal  A  =  -.01132093 

cpi  A  =  +.00689651* 

trl/k  =  -  .0021*9631* 
C63A  =  -.00905376 

t-U  A  -  -.00007979 
It]  A  -  +.00612325 
LAG  A  =  +» 01157250 
A  -  -.0028911*8 
U'fl  A  -  +.00721310 
[-L' 3  A  -  -.0000571*7 
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I 


The  observation  equations  ares* 


Aa2  "  100^2 

I5S4b2 

aV=  isb4®. 

Ab  "*»  JLfib 
0  100  c 

» 

A 

B 

C 

D 

-L 

(a»2) 

(Ab2) 

(Aao) 

(Abo) 

3 

+  .07.097260 

+.01.114795 

-.00.878273 

-.00.154878 

+.00002691 

10 

+.00.782872 

-.06.363547 

-.00.365933 

+.02.206392 

-.00001590 

17 

-.02.557324 

+.05.652257 

-.00.50540 3 

-.01.613615 

+ .00000031 

18 

-.05.323809 

-.00.403505 

+.01.749610 

-.00.437897 

-.00001133 

3 

+.01.522675 

-.06.930350 

-.00.388429 

+.00.989768 

+ .00000244 

10 

-.06.226825 

-.01.151650 

+.01.863052 

+.00.719163 

+ .00000361 

17 

+.05.419675 

+.02.858550 

-.01.226356 

+  .00.097364 

+.00002210 

18 

-.00.715525 

+.05.223450 

-.00.248267 

-.01.806294 

-.oooa28i4 

it 

The  decimal  point  in  the 

coefficients  A. 

..D  was  moved 

for  convenience 

in  the  numerical  computations. 


The  normal  equations  ares 


( *b'2 ) 


(AbV 


+11.6273018  +0.0301952 
+12.oo455o5 


-2.4368121 
-0.423?  $78 
+0.674/305 


+0.4123548  +.000027579431 
-2.9733715  +.000003680403 
+0.0342156  -.000003973645 
+0.9178321  + .000001339620 
+ .000000000181 


=  0 

e»  0 
=  0 

-  0 

-  0 


The  solution  of  the  normal  equation  gives? 


Aa  '  =  -.0000121886 
.0 

A.bo>  =  -.  0000062  844 
Aa2‘  =  -.0000046976 

Abg5  =  -.0000022811 

+.0000000000835 


AaQ  *  -.00121886 

AbQ  =  -.00062844 

Aag  =>  -.00046976 
Ab2  -  -.00022811 


And  with  the  formulas  (44)  and  (50)t 


Aax  =  -.00043917 

Abx  =  +.00001070 

Ac-j^  =  +.00006731 
ACg  =  +.00006198 
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The  plate  constants  are  nows 


a1  =  +2  <>22799080 

*  +2.57821*033 
c1  =>  +0.28205068 
a2  -  -2.71*21*9211* 
b2  -  +2.09U3M58 
c2  «  +0.22919672 
aQ  =  +0.00092915 
bQ  .  -1.20517253 


The  condition  equations  reduce  to: 
1)  +.00000090 
2)  +.0000011*6 


If  the  solution  is  repeated  to  eliminate  the  influence  of  second  order 
terms,  we  obtain  the  following  new  L  and  L*  terms: 


-L 

-L* 

3 

+.00000269 

+ .00000962 

10 

-.ooooii*i*5 

+ .00000827 

17 

+.00001571 

+  .000001*1*5 

18 

-.00000399 

-.00002331 

and  the  new  absolute  terms  in  the  normal  equations  ares 

[>(L)]  =  -.000000005030 
[B(L)]  -  -.000000022879 

[c(l)1  =  +.000000001219 

[D(L)]  *  +.000000006321* 

[L(L)1  =  +.000000000083 


The  corresponding  reduction  gives: 

Aa1  *  +.00000033 
^b]_  -.00000019 

AC1  =  +.00000001 

Aa2  =  +.00000016  . 

Ab2  =  +.00000007 

Ac_  =  -.00000001 

4.  .  +.000OO0U8  ■  ♦•000000000083U 

0 

Ab0  =  -.00000055 


And  the  final  plate  constants  ares 


&1  =  +2.22799078  *  a2  =  -2.7421*9198  aQ  -  +0.00092963 

bx  =  +2,57824013  b2  =  +2.09411*465  bQ  =  -1.20517308 

c±  =  +0.28205069  c2  =  +0.22919671 

The  condition  equations  reduce  now  to* 

1)  +.00000003 

2)  -.00000010 


The  computation  of  the  residuals  gives? 


P 

P 1 

V 

Vs 

3 

+2.709 

+9.600 

+1.68 

-2.32 

10 

-14.457 

+8.254 

+4.50 

+1.60 

17 

+15.726 

+4.444 

-1.54 

-3.80 

18 

-3.978 

-22.309 

-4.18 

+4.21 

834 

£w> 

831 

Mean  error  of  an  observed  1  or  l1  value  is: 


m  •  /H5  =  ~ 6*^ 

The  final  check  with  equation  (l4)  gives: 


I+v 

3  +.02135168 
io  -.05614050 
17  +.06031846 
18  -.00103618 


-.05773332 
+.00005760 
+.04015420 
+  .06381121 


< 

c 

V-V 

+.16900890 

8tfi  dec . 
+1 

+.04650154 

8thc  dec 
-1 

+.15713779 

0 

+.38332882 

-1 

+  .54637688 

0 

+.15537272 

-1 

+.48127492 

-1 

+.3^613273 

+1 

The  computation  of  the  orientation  elements  with  formula  (16)  gives: 


tan  K  = — ,00077137 
|  =  -0°  .2'  39.11" 


sin<=  -.00077137 
cos*  ®  +.99999970 
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7}±  =  ~„136372|60 
=  -.11095UU8 

cot  (A+<)  =  +lo2310li2UU 

(A+<)  =  219°  5*  I5.1i9'‘ 

sin  (A  +  *)  -  -.63050832 
cos  (A  +*)  =  776182^9 

A  =  219°  7 8  gk.59” | 

sin  A  «  -.6311068^ 
cos  A  =  -.77569591 

h*  =  +2o93l89h28 

X  '  »  +  .0000611*6 
c 

yc>  ■=  +.053099U6 
:  d'  =  +  .882855# 

sin  v  -  +.3U1076U2 

|  v  =  19°  5*6 8  32.69M 

v/2  «  9°  58’  16  ,W 

tan  v/2  =  +.17580936 

d  =*  ,30112121. 

y0‘  =  +,05293993 

yD'-ycJ  =  -.00015953 

Ax  =■  -.0000613^1 
Ay  =»  -.000159581 


compiled  results  of  the  three  adjustments  are: 


-3' O  Os  CM  eAM3  CO  CM 

o  d  g  e  e  o  o  o  *  6  ®  e 

iH  rH  CA-^  CM  |H  CA-O  CM  |H  CA-J 

0  +  JI  I+S+  a  +  0  + 


rH  C'-'O  CM 

«  0  o  e 

w  JrtJ 
+  +  0  ! 


UNCACO 
0*00 
pH  ?  rH 

+  H-  )  J 


t — "LA  U\  CM 
»  >  •  » 
pH*  -J  pH  -5 

+  +  a  8 


CA  O  C-CO  -CAO  C— CO  CAO  r— GO 


f- 

"LA 

rH 

■ 

« 

a 

Os 

f- 

CN 

CM 

CA 

Cr\ 

CM 

CM 

CNJ 

1 

a 

1 

ca 

CA 

Os 

e 

• 

o 

H 

CA 

CM 

ca 

CA 

<A 

SO 

M3 

VO 

U\ 

UN 

UN 

OS 

OS 

os 

pH 

pH 

rH 

C- 

M3 

M3 

e 

• 

e 

ca 

CA 

-O 

US 

UN 

>  c — 

c — 

f- 

Os 

Os 

Os 

£ 

rH 

CM 

& 

CM 

Os 

o 

SO 

US 

o 

3 

o 

pH 

8 

pH 

8 

o 

o 

o 

M 

0 

a 

a 

1 

a 

rH 

rH 

rH 

CM 

CVJ 

CM 

rH 

rH 

rH 

rH 

3 

o 

o 

O 

CA 

CA 

e 

Q 

• 

results  see  remark  on  page  73 


A  comparison  between  the  results  obtained  from  the  unique  3  stars 
solution  and  from  the  overdetermined  i*  stars  solution  shows  ’  ■/ 
an  important  fact.  The  square  sum  of  the  residuals  =  83. U  and  conse¬ 

quently  the  mean  error  of  an  individual  observed  plate  coordinate 


=  +  6,!?p>  or  expressed  in  angular  terms  with  a  plate  dis¬ 
tance  of  d~300  mm,  m  =  +  U".  We  should  now  expect  that  the  individual 
parameters  should  differ  only  by  a  comparably  snail  amount  between  the 
two  solutions.  However,  between  the  unique  and  overdetermined  solutions 
there  are  differences  of  as  much  as  £00**  for  the  individual  parameters. 

This  shows  that  the  individual  parameters  are  determined  with  relatively 
large  mean  errors.  Only  their  combined  effect  produces  a  high  precision 
result  in  representing  an  individual  ray  in  space. 

Consequently,  some  of  the  parameters  are  able  to  compensate,  at 
least  partly,  for  errors  on  the  other  parameters.  This  opens  the  way,e.g., 
to  determine  in  explicit  form  the  tolerances  for  a  camera  design  and  allows, 
on  the  other  hand,  a  shortening  of  the  numerical  computations  by  assuming 
some  of  the  unknowns  as  constants.  It  is  obvious  that  an  error  in  A  x  or 
Ay  will  be  partly  compensated  for  by  changes  in  azimuth  or  tilt, respectively. 
The  narrower  the  bundle  of  rays  the  more  effective  is  this  compensation. 

By  such  a  measurement  the  number  of  unknowns  will  be  reduced  to  four  and 
therefore  the  final  number  of  normal  equations  after  thS  elimination  of 
A  A  will  be  only  three. 

VII  SUMMAR? 


A  method  of  high-precision  spatial  triangulation  is  described.  The 
relative  angular  accuracy  of  the  method  is  better  than  1:100000.  The 
mathematical  analysis  of  the  plate  orientation  as  well  as  the  triangula¬ 
tion  by  spatial  resection  is  based  on  rigorous  geometry.  A  rigorous  least 
squares  adjustment  determines  the  most  probable  plate  orientation  and  de¬ 
livers  the  most  probable  corrections  to  the  measured  plate  coordinates, 
thus  making  it  possible  to  determine  the  mean  error  of  ;the  orientation 
elements  and  of  the  target  point  directions.  The  relative  angular  mean 
error  of  a  final  spatial  direction  to  any  recorded  target  point  is  1  to  2 
seconds  of  arc.  The  triangulation  is  based  on  the  principle  of  photogram- 
metry  by  intersection.  Hence,  each  point  is  determined  by  combining  the 
measured  results  of  individual  camera  stations .  The  mean  error  of  the 
final  coordinates  depends  on  the  geometry  of  the  individual  tri angulation 
case. 


i 
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The  method  is  useful  for  the  determination  of  the  spatial  position 
of  single  points ,  e.  g„,  for  the  purpose  of  calibration  of  other  trajectory 
measuring  methods,  as  well  as  for  the  determination  of  complete  trajectories. 
The  outlined  method  is  free  of  systematic  errors  if  the  conditions  discussed 
in  Chapter  V  are  sufficiently  satisfied.  It  must  be  understood  that  maxi¬ 
mum  accuracy  with  the  method  can  be  obtained  only  if  during  the  time  inter¬ 
val  between  the  exposure  of  the  control  points  for  plate  orientation  and 
the  exposures  of  the  target  points,  the  interior  and  exterior  orientations  of 
the  cameras  can  be  considered  as  unchanged* 


H.  Schmid 
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